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Chapter 1 


Introduction 


1.1. Assets 


Assets represent value of ownership that can be converted into cash. There are two kinds of assets: 
tangible and intangible. Commodity, foreign currency, house, building, equipments are tangible, 
while copyright, trademarks, patterns, computer programs and financial assets are intangible. 


The financial assets include bank deposits, debt instrument, stocks and derivatives. Debt instru- 
ments are issued by anyone who borrows money — firms, governments, and households. They in- 
clude corporate bonds, government bonds, residential, commercial mortgages and consumer loans. 
These debt instruments are also called fixed-income instruments because they promise to pay fixed 
sums of cash in the future. 


The stocks, shares and equities are all words used to describe what is essentially the same thing. 
They are the claim of the ownership of a firm. When someone buys a share, he is buying ownership 
of part of a company and becomes a shareholder in that company. 


The prices of assets are settled through trading. Many assets are traded in markets. There are 
commodity markets as well as financial markets. Financial markets include bond markets, stock 
markets, currency markets, financial derivatives markets, etc. The financial market provides a link 
between saving and investment. Savers can earn high returns from their saving and borrowers can 
execute their investment plans to earn future profits. 


The derivatives are contracts derived from some underlying assets. Usually, the prices of assets 
fluctuate time by time. In order to reduce the risk of price fluctuation, a corresponding contract 
is introduced to make such uncertainty more certain. For instance, suppose today’s price of corn 
is US$3 per bushel. You want to buy 5,000 bushels of corn for delivery two months later. You 
can sign an agreement with someone who is willing to sell you this amount at such price at such 
future time. This agreement is called a forward contract. It certainly costs you some money but 
make your price uncertainty more certain. We call that such forward contact hedges the risk of price 
fluctuation. Forward contract is one kind of financial derivative. The corn is the underlying asset. 
More financial derivatives will be introduced later. 
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1.2. Financial Derivatives 


1.2.1. Forwards contract 


A forward contract is an agreement which allows the holder of the contract to buy or sell a certain 
asset at or by a certain day at a certain price. Here, 


e the certain day—maturity or expiration date, 
e the certain price—delivery price, 
e the person who write the contract (has the asset) is called in short position, 


e the person who holds the contract is called in long position. 


Example 1. Suppose today’s price of corn is US$3 per bushel. You want to buy 5,000 bushels of 
corn for delivery two months later. You can sign an agreement with someone who is willing to sell 
you this amount at price, say US$3.05, at such future time. Then you make your uncertain risk of 
price fluctuation more certain. 


Example 2. (quoted from Chan) Suppose it is May 30 now and you need to pay your tuition 
to Oxford University £1,000 pound by September 1. Suppose you can earn HK$15,600 in the 
summer. The exchange rate now is £1.00 = HK$15.40. What should you do? 


e 1. Try your luck: If the exchange rate is less than HK$15.60 by August 31, you earn some 
extra cash besides the registration fee. If the exchange rate is higher than HK$15.60, you ... ? 


e 2. Buy a forward on British pounds: Look for a forward contract on British pounds that 
entitles you to use HK$15,600 to buy £1,000 on August 31. Then you have eliminated, or 
hedged, your risk. Thus, a forward contracts eliminate your risk. 


1.2.2 Futures (futures contracts) 


A futures contract is very similar to a forward contract. Futures contracts are usually traded through 
an exchange, or clearinghouse, which standardizes the terms of the contracts. The exchange helps 
to eliminate the risk of default of either party through a margin account. Each party has to pay 
an initial margin as deposit at the inception of the contract. The profit or loss from the futures 
position is calculated every day and the change in this value is paid from one party to the other. The 
maintenance margin is the minimum level below which the investor is required to deposit additional 
margin. 
The difference between forward contract and futures contract are 
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Forward Contract Future Contract 
nature customized contract | standarized contract 
Trading over the counter through exchanges 
liquidity less liquid highly liquid 
counter party risk | high negligible 
Settlement delivery (at the end) | closed out prior to maturity 
Margin no margin compulsorily needs to be paid by the parties 


Example = Mr. Chan takes a long (buy) position of one contract in corn (5,000 bushels) for March 
delivery at a price of US$3.682 per bushel at the Chicago Board of Trade (CBOT). See quotation 
atlhttp: //www.cbot.com/.CBOT|It requires maintenance margin of US$700 with an initial 
margin markup of 135%, i.e. the initial margin is US$945 which Mr. Chan and the seller each has 
to deposit into the broker’s account on the first day they enter the contract. The next day the price 
of this contract drops to US$3.652. This represents a loss of US$0.03 x 5,000 = US$150. The 
broker will take this amount from Mr. Chan’s margin account and deposit it to the seller’s margin 
account. It leaves Mr. Chan with a balance of US$795. The following day the price drops again to 
US$3.552. This represents an additional loss of US$500, which is again deducted from the margin 
account. As this point the margin account is US$295, which is below the maintenance level. The 
broker calls Mr. Chan and tells him that he must deposit at least US$405 in his margin account, 
or his position will be closed out, i.e. both sides agree to settle the contract at this point. Once the 
position is closed, Mr. Chan will not be able to earn back any money in the future even if the price 
rises above US$3.682. 


1.2.3 Options 


There are two kinds of options — call options and put options. A call (put) option is a contract 
between two parties, in which the holder has the right to buy (sell) and the writer has the obligation 
to sell (buy) an asset at certain time in the future at a certain price. The price is called the exercise 
price (or strike price). The holder is called in long position, while the writer is called in short 
position. The underlying assets of an option can be commodity, stocks, stock indices, foreign 
currencies, or future contracts. 

There are two kinds of exercise features: 


e European options : Options can only be exercised at the maturity date. 

e American options : Options can be exercised any time up to the maturity date. 
Notation 

e ¢ current time 

e T maturity date 


e S current asset price 
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e Sr asset price at time T 
e F strike price 
e c premium, the price of call option 


e r bank interest rate 


Examples An investor buys 100 European call options on XYZ stock with strike price $140. 
Suppose 


E = 140, 
S = 138, 
T = 2months, 
c = 5 (the price of one call option). 


If at time T, Sp > E, then he should exercise this option. The payoff is 100 x (Sp — FE) = 
100 x (146 — 140) = 600. The premium is 5 x 100 = 500. Hence, he earns $100. If Sp < T, then 
he should not exercise his call contracts. The payoff is 0. 

The payoff function for a call option is A = max{ Sr — E,0}. One needs to pay premium (c;) 
to buy the options. Thus the net profit from buying this call is 


A- cer F—-9) : 


Example Suppose 


todayist = 03/09/2016, 
expirationis 7) = 31/12/2016, 
the strike price F = 250 


for some stock. If S¢ = 270 at expiration, which is bigger than the strike price, we should exercise 
this call option, then buy the share for 250, and sell it in the market immediately for 270. The payoff 
A = 270 — 250 = 20. If S7 = 230, we should give up our option, and the payoff is 0. Suppose the 
share take 230 or 270 with equal probability. Then the expected profit is 

1 1 
Fg = et 
Ignoring the interest of bank, then a reasonable price for this call option should be 10. If Sp = 270, 
then the net profit= 20 — 10 = 10. This means that the profits is 100% (He paid 10 for the option). 
If Sp = 230 the loss is 10 for the premium. The loss is also 100%. On the other hand, if the 
investor had instead purchased the share for 250 at ¢, then the corresponding profit or loss at T’ is 
+20. Which is only +8% of the original investment. Thus, option is of high risk and with high 
return. 
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1.3 Payoff functions 


At the expiration day, the payoff of a future or an option is the follows. 


Payoff of futures 


e Payoff of a future in long position: At the expiration day, the price of the asset is Sp. The 
holder can buy the asset at price FE’. He has the obligation to buy it. Thus the payoff of the 
holder is Sp — E. 


e Payoff of a future in short position: at expiration, the writer also has the obligation to sell the 
asset at price E. So, he needs to use Sp to buy the asset, then sell to the holder at price EF. 
Thus his payoff is EF — Sr. 


(a) left (b) right 


Figure 1.1: Payoff of a future, long position (left) and short position (right) 


Payoff of call options 
e Long call: A = max{S — E,0} 
- If Sp > E, then the holder exercise the call at price F then sell at price Sp 
- If Sp < E, then he gives up the call option. 
e Short call: A = min{ E — Sy, 0} 


- If Sp > E, he has the obligation to sell the asset at price /, thus he needs to buy the 
asset at price Sv then sell it at price . He losses Sp — E. 


- If Sp < E, the holder will give up the call option. So the writer losses nothing. 
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call option (long): max{S,-K,0} call option (short):—max{S,-K,0}=min{K-S,,0} 


(a) left (b) right 


Figure 1.2: Payoff of a call, long position (left) and short position (right) 


Payoff of put options 
e Long put: A = max{ FE — Sr, 0} 


— If Sp < EF, the holder has the right to sell the asset at FE, then buy it back at Sip. Thus 
the payoff is & — Sr. 
- If Sp > E, the holder just gives up the option. 


e Short put: A = min{S7 — E,0} 


-— If Sp < E, he has the obligation to buy the asset at E. He then sell it at S7. Thus he 
losses EF — Sp. 


- If Sp > E, the holder will give up the put option, So the writer losses nothing. 


1.3.1. Premium of an option 


The person who hold an option has the right and the counter party has the obligation to fulfill the 
contract. Thus, a premium should be paid by the holder to the writer. Below is a portion of a call 
option copied from the Financial Times. 


the current timet = Feb3 
the expiration 7 = _ end of Feb, 
T—t & 10days 
S, = 2872 


E | 2650 2700 2750 2800 2850 2900 2950 3000 
c | 233 183 = 135 89 50 24 9 3 
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1.4 


1.5 


1.6 


put option (long):max{K-S,,0} put option (short):-max{K-S,,,0}=min{S,-K,0} 


(a) left (b) right 


Figure 1.3: Payoff of a put, long position (left) and short position (right) 


Other kinds of options 


Barrier option: The option only exists when the underlying asset price is in some prescribed 
value before expiry. 


Asian option: It is a contract giving the holder the right to buy or sell an asset for its average 
price over some prescribed period. 


Look-back option: The payoff depends not only on the asset price at expiry but also its 
maximum or minimum over some period price to expiry. For example, A = max{J— So, 0}, 
j= Maxo<7<T ST); 


Types of traders 


. Speculators (high risk, high rewards) 
. Hedgers (to make the outcomes more certain) 


. Arbitrageurs (Working on more than one markets, p12, p13, p14, Hull). 


Basic assumption 


Arbitrage opportunities cannot last for long. Only small arbitrage opportunities are observed in 
financial markets. Our arguments concerning future prices and option prices will be based on the 
assumption that “there is no arbitrage opportunities”. 
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Figure 1.4: The FT-SE index call option values versus exercise price. 


Chapter 2 


Asset Price Model 


2.1 Efficient market hypothesis 


The asset prices move randomly because of the following efficient market hypothesis: 


1. The past history is fully reflected in the present price, which does not hold any future infor- 
mation. This means the future price of the asset only depends on its current value and does 
not depends on its value one month ago, or one year ago. If this were not true, technical 
analysis could make above-average return by interpreting chart of the past history of the asset 
price. This contradicts to the hypothesis of no arbitrage opportunities. In fact, there is very 
little evidence that they are able to do so. 


2. Market responds immediately to any new information about an asset. 


2.2 The asset price model 


We shall introduce a discrete model and a continuous model. We will show that the continuous 
model is the continuous limit of the discrete model. 
2.2.1 The discrete asset price model 


The time is discrete in this model. The time sequence is nAt, n € N. Let us denote the asset price 
at time step n by S;,. We model the asset price by 


(2.1) 


Sn+1 u with probability p 
Sn 


~ | d_ with probability 1 — p. 


Here, 0 < d < 1 < u. The information we are looking for is the following transition probability 
P(S;, = S|So), the probability that the asset price is S at time step n with initial price So. We shall 
find this transition probability later. 
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2.2.2 The continuous asset price model 


Let us denote the asset price at time ¢ by S(t). The meaningful quantity for the change of an asset 
price is its relative change 

dS 

9” 
which is called the return. The change a can be decomposed into two parts: one is deterministic, 
the other is random. 


e Deterministic part: This can be modeled by 


dS 
—_— = t. 
Cie 


Here, js is a measure of the growth rate of the asset. We may think ju is a constant during the 
life of an option. 


e Random part: this part is a random change in response to external effects, such as unexpected 
news. It is modeled by a Brownian motion 


odz, 


the a is the order of fluctuations or the variance of the return and is called the volatility. The 
quantity dz is sampled from a normal distribution which we shall discuss below. 


The overall asset price model is then given by 


< = pdt + adz. (2.2) 


We shall look for the transition probability density function P(S(t) = S|S(0) = So). Or equiva- 
lently, the integral 


[rs = $|S(0) = So) dS 


is the probability that the asset price S(t) lies in (a, b) at time t and is So initially. 


2.3. The solution of the discrete asset price model 


Let us consider the discrete price model 


Snti — f wu with probability 1/2 
Si d_ with probability 1/2. 


A sequence of movements (.59, 51, ..., Si.) is called an n-step path. In such a path, it can consist of 
£ times up movements and n — @ times down movements, where 0 < & < n. The corresponding 
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values of S,, are Syu“d"—*. Since the probability of each movement is independent, the probability 
of an n-step path with ? up movements (and n — € down movements) is 


1 £ 1 n—l 7 1\” 

2) \2 “Kp 
There are (7) paths with @ up movements in n-step paths, we then obtain the transition probability 
of the asset price to be: 


(7) (4)” when S = Spu’d”—£, 


0 otherwise. G2) 


P(Sn = SIS) = { 


Such S,, is called a (discrete) random walk. It is hard to deal with the values Sou“d"—*. Instead, we 
take its logarithmic function. Let us denote In S,, by X;,. This random walk X,, obeyes the rule 


XY 4-X= Inu with probability 1/2 

nyt" ') Ind with probability 1/2. 

Let us define i ind , ind 
Ar=———, In(1 +r) = SSE. 


Then the rule can be written as 


a AG with probability 1/2 
Angi a= ihe) { —Ax with probability 1/2. 
The term In(1 +1) * r is called a drift term. It measures the growth of S,,. It can be absorpt into 
Xn. Indeed, if we define Z,, = In((1+ r)~"S,,), then Z,, satisfies 


Ag with probability 1/2 


Znt — 4n = { —Ax with probability 1/2. ea) 


The advantage of this formulation (using Z,, now) is that this random walk has equal increment Ax, 
the possible values of the random variable Z,, are mAz, m € Z. This is the random walk with 
binomial distribution. 


2.4 Binomial distribution 


Consider a particle moves randomly on an one dimensional lattice {mAz|m € Z} according to 
the rule (2.4), where Z,, denotes the location of this particle at time step n. Suppose the particle is 
located at 0 initially. Our goal is to find the transition probability 


Pl 4, =| 26 =0) 


explicitly and also find its properties. 
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This random walk is equivalent to the following binomial trials: flipping a coin n times. Suppose 
we have equal probability to get Head or Tail in each toss. We move the particle to the right or left 
adjacent grid point according to whether we get Head or Tail. Let X;, be the result of the kth 
experiment. Namely, 

= { 1 if we get Head 
—1_ if we get tail 


Then Z,,/Ax = yoy X;. The random walk Z,, is equivalent to the random walk Y,,, the number 
of Heads we get in 7 trials. Namely, 


PZ, = mel Zo= 0) =] PY, =2) 


with j 
m =—(n—£) = 2€—n, or equivalently ¢ = a(n +m). 


Notice that m is even (odd), when n is even (odd). There is a one-to-one correspondence between 
{Z\0<l<nbo{m| —-n<m<n,m+ nis even}. 


Thus, we can also use £ instead m to label our particle. There are CA paths to reach mAzx from 0 
with 2 = (n+ m)/2, we get 


0, if m + nis odd, 
(arenes) (5)”, if m+ 7 Is even. 


ros=0=()() 


£=0 


This probability distribution is called the binomial distribution. We denote them by pz, and py,, 
respectively. 


Or equivalently, 


It is clear that 


Properties of the binomial distribution We shall compute the moments of pz, or py,,. They are 
defined by 


ie ee a a) 
¢=0 m= on 


m + neven 


Since m = 2¢ — n, we can find the moments < m* >=< (2¢—n)* > by computing < @* >, 
which in turn can be computed through the help of the following moments generating function: 


cn 4mo EH) ()-CEY- 
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To compute moments, for instance 


G'(1) = ye - py, (0) =<k>. 
¢=0 


n 


G"(1) = 5) e(é- Dpy, (2) =< P > - <e>. 


l=1 
Hence, we obtain the mean 
n 
aL>=G'il) = = 
The second moment 
n(n + 1) 


<0 S=G" (+ <ts= 7 
We can also compute the mean and variance of pz,, as the follows. From m = 2¢ — n, we have 
moa) <= (oan =U): 

To compute the second moment < m? >, from m = 20 — n, we have 
<m>=4<0>—-4n<€>4n?% =n. 

The mean of this random walk is < m >= 0, while its variance is < (m— <m >)? >=n. 

Exercise 

1. Find the transition probability, mean and variance for the case 

ae ee { AG with probability v) 

—Az_ with probability 1 — p 


2.5 Normal distribution 


A random variable Z is said to be distributed as standard normal if its probability density is 


= 1 —2? [2 


pz(x) = ia 


We denote it as Z ~ N’(0, 1). A random variable X is said to be distributed as normal with mean ju 
and variance o° if its probability density function is 


oo (@—H)?/20? 


1 
PX\%) = 
o Vv 210? 
We denote it by X ~ N (ju, 07). For such .X, it is clearly that 


= 
Li= ad 


~ N(0,1). 
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For 

X— pL oe 1 2 ND 2 ane =D 
Pla< <b) = P(utoa < X < ptob) = ee es W180) de = | ae de 
( ‘ ptoa VIro2 a Jor 


= 


The moments of Z (i.e. E/Z™ 


]) can be computed through the helps of moment generating function 


e° s™ 7 mm 
G(s) := a -> a fe” a pz(e)de= y ra u|Z”™ |. 
m=0 
We have 
CO ok 
est ene? /2 = [2 l / —(a—s)?/2 _ af /2 $ 
G(s) = dz =e Tae. dz =e d DEF 


By comparing the two expansions above, we obtain 


2") = 0 if misodd 


Ch ifm = 2k. 


2.6 The Brownian motion 
2.6.1 The definition of a Brownian motion 
The definition of the Brownian motion (or the standard Wiener process) z(t) is the following: 
Definition 2.1. A time dependent function z(t), t € R is said to be a Brown motion if 
(a) Vt, z(t) is a random variable. 
(b) z(t) is continuous in t. 
(c) For any u > 0, s > 0, the increment z(t + s) — z(t), z(t) — z(t — uw) are independent. 
(d) Vs > 0, the increment 2,445 — 2, is normally distributed with mean zero and variance s, i.e. 


Zt+5 — 2 ~ N(0,8). 


2.6.2. The Brownian motion as a limit of a random walk 


We may realize the Brownian motion as the limit of the standard random walk. Let us study Brown- 
ian motion in [0, ¢]. First, we partition the time interval [0, t] into n subintervals evenly. Let us imag- 
ine a particle moves randomly on lattice points {mAz|m € Z} at discrete time kAt, k = 0,...,n 
We choose Ax = At. The motion of the particle is: 
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Zk = Zh-1 + X,Ax 

where 

x,- 1 with probability 1/2 

k~\ =1 with probability 1/2 

The location Z,, at time step n obeys the binomial distribution pz, (m) = P(Zn, = mAa|Zo = 0). 
We can connect discrete path (Zp, Z1, ..., Zn) by linear function and form a piecewise linear path 
Z"(s),0< 5s < t with Z"(kAt) = Z,. As n — ov, we expect that these kinds of paths tend to 
a class of zig-zag paths z(-), called Brownian motions. They are time-dependent random variable. 
Their properties are: 


e 2(0) =0; 

e z(-) is continuous; 

e Any non-overlapping increments z(t4) — z(tg) and z(t2) — z(t1) (ti < te < tz < ta) are 
independent. 


Let us denote the collection of these “zig-zag” paths by Q, the sample space of the Brownian mo- 
tions. 

Next, we find the probability distribution of Brownian motions. The probability mass function 
of the random walk is P(Z, = mAz|Zp = 0}. Since n + m is even, the corresponding probability 
density function on the real line is P?, :-= P(Z, = mAz|Zp = 0)/(2Az). We claim that 


Proposition 2.1. The probability density function of the random walk has the limit 


pry e~2?/(2t) 


1 
V 2at 
asn — co with mAx > x and nAt = t. 


This is the probability density of the standard normal distribution \’(0, t). Thus, the 4th property 
of the Brownian motion is 


e The probability distribution of z(t) is z(t) ~ N’(0, t). 
We prove this proposition by using the Stirling formula: 
nl x V2nn™tae™, 


Recall that the probability 


a ee a) (3) 


Using the Stirling formula, we have forn, m,n —m >> 1, 


(am'+m) (3) 


Al 
V2rn"t se” 


1 n 
(5) V2n(4(n + m))20+™) +3 e-(n-+m)/2 I (4(n — m)) 2 ™)+3 e—(n—m)/2 
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We shall use the formula seats 
(1+ =) >e", asn—>oo. 
n 


We may treat m as a real number because the functions above are smooth functions. Now, we fix 
x,t, define m = x/Ax = VCn, where C = x? /t. We take n — oo. Then 


(a) 2) se) ae ae 
+=) pr ( + | — ee? 
(2 _ oo“ = (: — 4 “ eT 


Oey G-Sy "OR sys 


Thus, we obtain 


We also have 


2Ar V2rnV/ At ~ V2nt 


Summarizing above, we obtain 


1 n ry" 1 2 
Pia. = —2@ /(2t)_ 
me De Cees (5) > Vint. 


Remarks. 
1. We can also show this limiting process as a corollary of the central limit theorem. Recall that 


n 
Zn = > X,Az, 
k=1 


where X;, ~ X are iid and 


1 with probability 1/2 
—1_ with probability 1/2 


be 
II 
——s 
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Since Ax = VAt = \/t/n, we get 


1 n 
Zn =—=Y VtXzp. 
a7? > 


The mean < /tX >= 0, the variance < (JEX - >= t. Thus, by the central limit theorem, 
the random variable Fi yp VtX» converges in distribution to a normal (0, t): 


1 n 
ae —+ N(0,t). 


2. We can consider more general random walk 


n 
Zn = > XpoAz, 
k=1 


where X;, ~ X are i.i.d. The parameter o, called volatility, measures the variation of Z in 
each movement. In this case, 


_i¢ 2 
Zn = nee —+ N (0,07). 


2.6.3 Properties of the Brownian motion 


By (d) of the definition 
1 a 
P (z(t) = x|z(0) = 0) = — =e 2. 
(2(t) = ale(0) = 0) = = 
Let us consider the Brownian motion starting from y at time s and reaches z at later time ¢. The 
probability density is denoted by p(s, y, x, t), called the transition probability density of the stochas- 


tic process z(t): 


p(s, y,t, 2) = P(z(t) = al2(s) = y) = P(2(t) = x — y|z(s) = 0) 
= P(z(t — s) = 2 —y|z(0) =0) = =o 


Let us simplify the notation p(0,0,t, x) by p(t, x). This transition probability density is the proba- 
bility density function for Brownian motions z(t) starting from 0 at time 0. Any function f(z), its 
expectation is defined to be 


ale) = f ™ f(2)p(t,2) de. 


This transition probability function has the following properties. 
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e7@/(2t) de — 0, 


ca dl 
a V 2rt 


a 1 
/ a2? ——e*/ 24) dy = t. 


co | WV2ant 


3. Let us consider an infinitesimal change dz := z(t + dt) — z(t). Here, dt is an infinitesimal 
time increment. Then we have 


(dz)? = dt with probability 1. (2.5) 


This formula is very important for the stochastic calculus below. It is interpreted by integrat- 
ing the above equation in t: 
t 
| (dz)? 4. 
0 


The left-hand side is defined to be the limit of the Riemann sum )77_, (2(te) — 2(tk—1))?. 


Proposition 2.2. Let z(-) be the Brownian motion. Then 


n 


Jim S"(2(te) —z(tp—1))* =t with probability 1 
k=1 


asn — oo. Here ty := kt/n. 


Proof. 1. We shall apply the strong law of large numbers. Let us consider the random variables: 
Y= n(z(te) — z(th—1))?. We see that 


n 


3 (e(te) ~ 2Ctea))? = =o Mi. 
k=1 


k=1 


The random variables Y;, are i.i.d. Y; ~ Y, where 
Y = n(z(t/n) — z(0))?. 


2. The mean of Y is 


2? /(2t/9) dep = nit/n) St. 


1 
E|Y] = [ea 
/2n(t/n)? 
The variance of Y is also finite (check). That is 


E|(Y — t)?] < c, (2.6) 


where 
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3. By the strong law of large numbers 


1 n 
— y Y; — t with probability 1. 
n 

k=1 


Exercise 


1. Find E[z(t)”]. 


2. Check (2.6). 


2.7. Itd’s Lemma 


In this section, we shall study differential equations which consist of deterministic part: « = b(x), 
and stochastic part 02(t). Here, z(t) is the Brownian motion. We call such an equation a stochastic 
differential equation (s.d.e.) and expressed as 


dx(t) = b(x(t))dt + o(a(t))dz(t). (2.7) 
An important lemma for finding their solution is the following It6’s lemma. 


Lemma 2.1 (It6). Suppose x(t) satisfies the stochastic differential equation (2.7). and f (x,t) isa 
smooth function. Then f (x(t), t) satisfies the following stochastic differential equation: 


1 
df = (4 Abit 50° fox dt + of, dz. (2.8) 
Proof. According to the Taylor expansion, 


1 1 
df a fr dt + fx dx aa af (dt)? + ft dx dt + gfe (dx)? Se 
Plug into this equation. The term 
(dx)? = b? (dt)? + 2bodt - dz + 07(dz)?. 


In the Taylor expansion of df, the terms (dt)°, dt- dz are relative unimportant as comparing with the 
dt term and (dz)? term. Using and noting (dz)? = dt with probability 1, we obtain (2.8). 
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Example 1 The stochastic process 
x(t) := x9 + at + oz(t) 
solves the s.d.e. 
dx = adt + adz, 


where a and o are constants. By letting y = x — at, then y is a function of x and ¢. From Ité’s 
lemma, y(a(t),t) is also a stochastic process and y satisfies dy = odz. Since o is a constant, we 
then have y(t) = yo + oz(t). The transition probability density function for y is 


1 
V2r0*t 


Or equivalently, the transition probability density function for x is 


eo (Y— Yo)? /207t_ 


P(y(t) = yly(0) = yo) = 


1 
P(a(t) = a|x(0) = x0) eS en (e—at—a0)? /207t 


V2n0%t 
Example 2. Let y(t) = yo + 2(t)?. Then y(t) solves 
dy = dt + 2z(t)dz(t). 
To show this, we apply It6’s lemma. Let y = f(a, t) = 2”, and 2 = z. Then 
dy(t) = frdx + ; ferdt = 2z(t)dz(t) + dt. 


This means 
d(z*) = Qzdz + dt. 


Example 3. Let dx = dz and y = Inz. Then f, = 1/z and fy = —1/x?. Thus, we have 
1 1 
In terms of Brownian motion, it can be written as 
dlnz = — — — . 
z 
Thus, In z solves the above s.d.e. 
Example 4. Let y = e°”, o is aconstant. We choose dx = odz and y = e*. Then 
1 1 
de” = e* dx + 50 (dex)? = e"odz+ x0 at. 


Thus 1 
de°’* = e"*odz + 50 dt. 
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2.8 Geometric Brownian Motion: the solution of the continuous asset 
price model 


In this section, we want to find the transition probability density function for the continuous asset 
price model: 


dS = pS dt + oS dz. (2.9) 


with initial data S(0) = So. We apply It6’s lemma with x = f(S) = In S. Then, fs = 1/5 and 
fss = —1/S?. By Ité’s lemma, z satisfies the s.d.e. 


o2 


1 
dx = fgdS + 5fss(asy° = (u = =) dt +o dz. 


and x(0) = xo := In Sp. From Example | of the previous section, we obtain 


o2 
o(t) = 40 + (u— Ft + o2(0) 


Since S' = e”, we obtain 
2 
S(t) = Soexp ( _ se + a2(t)) ‘ (2.10) 


The probability density of x(t) is 


1 a2 \4)2 2 
Pre) = 2|2(0) S79) = en (e—to- (uF Yt)? /207t 
(x(t) = x|x(0) = 20) iar 
That is, 
a 
x(t) ~ N (ato + (u- at o*t) 
From 


P(S(t) = S|S(0) = So)dS = P(ax(t) = z|z(0) = x0) de 
=P(r(t) = a\e(0) = a29)dS/S 


we obtain that the transition probability density function for S(t) is 


1 hie 8 Gye. = 
P(S(t) = S|S(0) = So) = Vanos ggg hoa ue (2.11) 


This is called the log-normal distribution. The stochastic process x(t) obeys the normal distribution, 
while S(t) = e* is called to obey the geometric Brownian motion. 
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Remark. Alternatively, we may also use probability distribution function to derive the probability 
density function as below. The probability distribution function of S(t) is 


S 
Fs(S) = i? P(S(#) = Si|S(0) = 5) 4S, 


On the other hand, 


We notice that 


Thus, 


Figure 2.1: The log-normal distribution. 


Properties of the geometric normal distribution S(t) Let us denote the transition probability 
density P(S(t) = s|S(0) = So) by psi)(s) for short. That is 


1 (in £ (ut)? /207t 
Ps(t)(s) = a paere : 4 


We compute its mean and variance below. 
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Proposition 2.3. The mean and variance of the geometric Brownian motion are 


(a) E[S(t)] = J°>, spsq(s) ds = Soe, 


(b) Var|S(t)] = S2e2#[e?"t — 1). 


Proof. (a) 
ia etn 8 uF)? /2074 
E[ S(t =i 8 sas = | € So 2 ds 
[ (t)] H Pst) ( ) f a Ora 2k 
"2 1 a2 \4)2 jy 2 
=| er e7 (2-0 (uF )t) /20 t dex 
-~co V 2ra*t 
ee 1 22 4\2 jo 2 
-| etttotut.—(@+ Ft) [20 t dx 
—co V2r07t 
a u ~ («+ t)? /20%t 
= Soe" ‘| a PP ae 
—coo V2r0t 
v2 1 22 4\2 Jo 2 
= Sel" | e @- Ft) [20 t dr 
—oo Qro2t 
= Soet. 
(b) 


B[S(t)"] = [ s*psit)(s) ds = is oH FF 207 
0 0 V2n02t 
—_ is 1 o2® p—(#—20-(u- Ft)? /207t ie 


~co V2702t 
i : 2(a+20-+t) (a+ %-t)?/20%t 
= € e€ 2 dx 
—co V2702t 
22 [~ 1 2n—(x-+2~t)2/2 2t 
= Soe e 2 ile 
—oo Qro2t 


2 1 302 4)2 jo 2,1 2 

_ Sel" | e  @- Fh) /20*t+o t dx 
co V27r0"t 

_ @2 2ut+o7t 

= Soe ; 


Var[S(t)] = E[S(t)2] — E[S(t)]? = S2e2"[e"t — 1]. 


2.9 Calibrating geometric Brownian motion 


How do we obtain yz and o of an asset? We estimate them from the historical data 5;, at discrete 


time t; := iAt, i = 0,...,n. The procedure is as below. First, let x(t) = In S(t). Let Ax; := 
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x(ti41) — x(t;). The variables Ax; are independent random variables with identical distribution: 
a2 
Ag ~ At~ N((u- At 7 At), 


The historical data we collect are S; which is a sample of S(t;). We define x; = In S;. We estimate 
the mean and variance of Ax by using these historical data «;: 


6 cae a 
amy -_ Tj4+1— Mj) = ™M, 


i=1 


The fz and G are the estimators of wand o. 


'The variance estimator has n — 1 in denominator instead of n. This is called unbiased estimator because the expec- 
tation of this estimator is the correct variance. This is called Bessel’s correction. 


Chapter 3 


Black-Scholes Analysis 


3.1 The hypothesis of no-arbitrage-opportunities 


The option pricing theory was introduced by Black and Scholes. The fundamental hypothesis of 
their analysis is that there is no arbitrage opportunities in financial markets”’. 
For simplicity, we shall also assume 


1. There exists a risk-free investment that gives a guaranteed return with interest rate r. ( e.g. 
government bond, bank deposit.) 


2. Borrowing or lending at such riskless interest rate is always possible. 
3. There is no transaction costs. 

4. All trading profits are subject to the same tax rate. 

We will use the following notations: 


current asset price 

exercise price 

expiry time 

current time 

growth rate of an asset 
volatility of an asset 

asset price at T 

risk-free interest rate 

value of European call option 
value of American call option 
value of European put option 
value of American put option 
the payoff function 


PyS QT HIE THER 
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Remark. 


1. In assumption 1, if the annual interest rate is r and the compound interest is counted daily, 
then the principal plus the compound interest in one year are 


AGA r 365 Ae" 
(tae) SaaS 


2. Assumption | implies that an amount of money A at time t has value Ae"? at time T if it 
is deposited into bank. This is called the future value of A. Similarly, a strike price F at time 
T should be deducted to Ee~"?—®) at time t. This is called its present value. 


3. The option price c, p, C and P are functions of (.S;, ¢). Usually, we only write c, which means 
the current value c(.S;,t). We use the same convention for S, an abbreviation of the current 
value S;. 


3.2 Basic properties of option prices 


3.2.1 The relation between payoff and options 


1. The payoff of a contract is the return from the deal. For European option, it can only be 
exercised on expiry date. Therefore, its payoff A is only meaningful at the expiry date T. 
However, for American option, the option can be exercised at anytime between ¢ and the 
expiration date JT’. It is therefore meaningful to define the payoff function at t. For a person 
who longs an American call option, his payoff is A(t) = max(.S; — E, 0), while for American 
put option, A(t) = max(E — S;,0). 


2. c(Sr,T) = A(T) = max(Sr — E,0). 

Otherwise, there is a chance of arbitrage. For instance, if c(S7,7T) < A(T), then we can buy 
a call on price c, exercise it immediately. If Sp > FE, then A = Sp — E > Oandc < A by 
our assumption. Hence we have an immediate net profit Sp — E — c > 0. This contradicts 
to our hypothesis. If Sr < E, then A(T’) = 0 andc < A(T) = 0 leading to c has negative 
value. This means that anyone who can buy this call option with negative amount of money. 
This again contradicts to our hypothesis. 

If c(Sr,T) > A(T), we can short a call and earn c(S7,7T). If the person who buy the call 
does not claim, then we have net profit c. If he does exercise his call, then we can buy an asset 
from the market on price Sr and sell to that person with price E. The cost to us is Sp — E. 
By doing so, the net profit we get is c(S7, T) — (Sr — E) > 0. Again, this is a contradiction. 


3. p(S7,T) = A(T) = max(E£ — Sr,0). 
If p(Sr,T) < A(t), then an arbitrageur longs a put at T, exercises it right away, and get net 
profit A(T’) — p(Sr,T). If p(Sr,T) > A(T), then an arbitrageur shorts a put. If the person 
who buys the put exercise it, then the arbitrageur earns P(S, 7) — A(T); if he gives up this 
put, then the arbitrageur earns p(S7, T). 
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4, We show C(S;,t) = A(t) := max(S; — E,0). 
If C(S;,t) > A(t), then an arbitrageur can short the call and earn C(S;, t). If the person who 
exercises it right away, the arbitrageur need to pay A(t). He still has net profit C(.S;, t) — A(t). 
If C(S;,t) < A(t). Then an arbitrageur buy C and exercises it right away and get net profit 


A(t) — C(S:,t). 


5. P(S;,t) = A(t) := max(E — S;,0). This can be shown similarly. 


3.2.2 European options 


Theorem 3.1. The following statements hold for European options 


max{S — Ee"F-®) 0} <e<3$ (3.1) 
max{Ee~"\7-9 — §,0} <p < Ee"? 4) (3.2) 

and the put-call parity 
p+S=ctEert (3.3) 


To show these, we need the following definition and lemmae. 
Definition 3.2. A portfolio is a collection of investments. 
For instance, a portfolio J = c — AS means that we long a call and short A amount of an asset 


o. 


Remark. The value of a portfolio depends on time. Suppose a portfolio J has A amount of cash 
at time t, its value is Ae"(T—*) at time T. On the other hand, a portfolio involves EF’ amount of cash, 
its present value is Ee~"7—), We should make such deduction, otherwise there is an arbitrage 
opportunity. 


Lemma 3.2. Suppose I(t) and J(t) are two portfolios containing no American options. Then under 
the hypothesis of no-arbitrage-opportunities, we can conclude that 


I(T) < J(T) > 1(t)< J), Ve<T, 


Proof. Suppose the conclusion is false, i.e., there exists a time t < T such that [(t) > J(t). An 
arbitrageur can buy (long) J(t) and short (t) and immediately gain a profit J(t) — J(t). Its value 
at time T is (I(t) — J(t))e"7—. Since I and J containing no American options, nothing can 
be exercised before T. At time T,, since J(T) < J(T), he can use J(T’) (what he has) to cover 
I(T’) (what he shorts) and gains another profit J(T) — I(T’). This contradicts to the hypothesis of 
no-arbitrage-opportunities. 

The equality case can be proven similarly. 
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Proof of Theorem[3.1) 
1. Let! =cand J = S. At T, we have 
I(T) = cr = max{Sr — E,0} < max{Sp,0} = Sp = J(T). 


Hence, I(t) < J(t) holds for all t < T. 
Remark. The equality holds when & = 0. In this case c = S 


2. Consider J = cand J = 0. At time T, J(T) = A(T) > 0 = J(T), hence I(t) > 0 for all 
t < T. Similarly, we also get p > 0. 


3. Consider J = c+ Ee~"7—4 and J = S. At time T, 
I(T) = max{Sr — E,0} + E = max{S7, FE} > Sr = J(T). 
This implies I(t) > J(t). 
4. Let I = pand J = Ee"? At time T, 
I(T) = max{E — S7,0} < E=J(T). 
Hence, p(t) < Ee~"7—). The equality holds only when Sy = 0. 
5. Consider J = p+ S and J = Ee~"7-), At time T, 
I(T) = max(£ — S7,0)+ Sr = max{S7, E} > E=J(T). 
Hence, I(t) > J(t). 


6. For put-call parity, we consider I = c+ Ee“ and J = p+ S. At time T, 


I(T) 
J(T) 


c+ EF = max{Sr — £,0}+ EF = max{S7, EF}, 
pt+S =max{E — Sr,0} + Sr = max{E, Sr} 


Hence, I(t) = J(t). 


3.2.3. Basic properties of American options 
Theorem 3.2. For American options, we have 
(i) The optimal exercise time for American call option is T, and we have C' = c. 


(ii) The optimal exercise time for American put option is as earlier as possible, and we have 
P(t) > p(t). 


(iii) The put-call parity for American option is 


C+ Bet) 29 Pe Oa. (3.4) 
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To prove these properties, we need the following lemma. 


Lemma 3.3. Let I or J be two portfolios that contain American options. Suppose I(r) < J(r) at 
some T <T. Then I(t) < J(t), for allt <r. 


Proof. Suppose I(t) > J(t) at some t < 7. An arbitrageur can long J(t) and short I(t) at time t to 
make profit J(t) — J(t) immediately. At later time 7, he can use J(7) to cover [(7) with additional 
profit J(7) — I(r), in case the person who owns I exercises his American option. 


Remark. The equality also holds , that is, if J(7) = J(r) for 7 < T, then I(t) = J(t) fort <7. 


Proof of Theorem|3.2} 


1. Firstly, we show C; > c; for all t < T. If not, then c(r) > C(r) for some time 7 < T, we 
can buy C and sell c at time 7 to make a profit c(7) — C(7). The right of C' is even more than 
that of c. This is an arbitrage opportunity which is a contradiction. 


Secondly, we show c; > C;. Consider two portfolios I(t) = C; + Ee~"\7~ and J(t) = St. 
The portfolio J is an American call plus an amount of cash Ee~"(7—), Suppose we exercise 
C at some time t < T, then I(r) = S; — E + Ee~"'—7) and J(r) = S;. This implies 
I(r) < J(r). By our lemma, I(t) < J(t) for all t < 7. Since 7 < T arbitrary, we conclude 
I(t) < J(t) for all t < T. Combine this inequality with the inequality 


ee be 4 Ss, 


in section 3.2,we conclude c¢ = C; for all t < T. Further, early exercise results in C(7) + 
Ee"-7) = §(r) — E+ Ee“"-7) < S(r). But S(r) < c, + Ee7"(7-7). Thus, early 
exercise leads to decrease the value of C(r), if 7 < T. Hence, the optimal exercise time for 
American option is 7’. 


2. We show p(t) < P(t). Suppose p(t) > P(t). Then we can make an immediate profit by 
selling p and buying P. We earn p — P and gain more right. This is a contradiction. 


Next, we show that we should exercise American puts as early as possible (the first time 
E > S,). We consider a portfolio I(t) = P, + S;. This is called a covered put. That is, 
we use S' to cover P when we exercise P. When we exercise P at T < T’, the payoff of I 
is E—S,+ 5, = E. By the time T’, its value is Ee'(?-7), Thus, in order to maximize the 
profit, we exercise time should be as earlier as possible. 


3. The inequality C + Ee~"\7—) < §$ + P follows from the put-call parity and the facts 
that c = Cand P > p. To show the inequality S + P < C+ E, we consider two portfolios: 
I(t) = c(t) + E and J(t) = P(t) + S;. Suppose P is not exercised before T. Then 
J(T) = max(£ — S7,0) + Sr = max(E, Sr), and 


I(T) =c(T) + Ee"? -® = max(E, Sr) + E(e"?-9 — 1) > I(T). 
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Suppose P is exercised at some time 7, with 7 < IT’. Then we sell S; at price E. Thus 
J(7) = E. On the other hand, 


I(r) =c(r) + Ee") > E= Jr). 


From lemma, we have I(t) > J(t). That is, c(t) + E > P(t) + S;. Since c(t) = C(t), we 
also have C(t) + E > P(t) + S;. 


—r(T-t) 


Remark. P — pis called the time value of a put. The maximal time value is F — Fe : 


Examples. 


1. Suppose S(t) = 31, EF = 30, r = 10%, T — t = 0.25 year, c = 3, p = 2.25. Consider two 


portfolios: 


I = c+ Ee? = 3430 x €701*95 = 39.96, 
J = p+tS=2.254+31 = 33.25. 


We find J(t) > I(t). It violates the put-call parity. This means that there is an arbitrage 
opportunity. 

Strategy : long the security in portfolio J and short the security in portfolio J. This results a 
cashflow: —3+ 2.25+31 = 30.25. Put this cash into a bank. We will get 30.25 x e9:!*9.25 = 
31.02 at time JT. Suppose at time 7, S7 > E, we can exercise c, also we should buy a 
share for F to close our short position of the stock. Suppose S7 < E, the put option will be 
exercised. This means that we need to buy the share for FE to close our short position. In both 
cases, we need to buy a share for EF to close the short position. Thus, the net profit is 


31.02 — 30 = 1.02. 


. Consider the same situation but c = 3 and p = 1. In this case 


T = et bet?) 33395 


and we see that J is cheaper. 

Strategy: We long J and short J. To long J, we need an initial investment 31 + 1, to short c, 
we gain 3. Thus, the net investment is 31 + 1 — 3 = 29 initially. We can finance it from the 
bank, and we need to pay 29 x e9-1*0.25 _ 99 73 to the bank at time T. Now, at T, we must 
have that either c or p be exercised. If Sp > EF, then cis exercised. We need to sell the share 
for £ to close our short position for c. If Sp < E, we exercise p. That is, we sell the share 
for &. In both cases, we sell the share for F. Thus, the net profit is 30 — 29.73 = 0.27. 
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3.2.4 Dividend Case 


Many stocks pay out dividends. These are payments to shareholders out of the profits made by the 
company. Since the company’s wealth does not change after paying the dividends, the stock price, 
the strike prices fall as the dividends being paid. If a company declared a cash dividend, the strike 
price for options was reduced on the ex-dividend day by the amount of the dividend. 

The inequalities below are all at the present time t. Thus, let D be the present value of the 
dividend, no matter when they are paid. Certainly its value at time Tis De’). Thus, a company 
pays dividend before T, its asset at time T becomes Sp + De"7—"), 


Theorem 3.3. Suppose a dividend will be paid during the life of an option. Let D denotes its present 
value. Then we have for European option 


Sap= bet eens (3.5) 
=SeD be YM <p < bet, (3.6) 
and the put-call parity: 
c+ Ee?) =p4+S—D. (3.7) 
For the American options, we have (i) 
S=DHPeO-Pe sabe (3.8) 


provided the dividend is paid before exercising the put option, or (ii) 
SHE SOS Pe espe! (3.9) 
if the put is exercised before the dividend being paid. 


Proof. 1. Proof of c> S — Ee~"\?-4) — D. 
We consider two portfolios: 


I = c+D+Ee7?-9, 
S. 


Then at time 7’, 


max{Sp — E,0}+ D+ E =max{Sp, E} + De" 
Sr t+ Det(F-*), 


NY 
—— 
4 
4 
| 


WD) 


Hence I(T) > J(T). This yields I(t) > J(t) for all t < T. This proves 
e=SaD=Fe tT, 


In other word, c is reduced by an amount D. 


32 


CHAPTER 3. BLACK-SCHOLES ANALYSIS 


2. Proof of p > D + Ee~"(7—-9) — S. That is, p increases by an amount D. 


We consider 
Lp 8. 
J=Ee*T-) + D. 


Attime T, I(T) = p(T) = max(E — Srp,0)+Sp+De"?, while J(T) = E+ De", 
Thus, I(T) > J(T), we then obtain p+ S > Ee~"(7-9 + D. 


. For the put-call parity, we consider 


I e+D+4+ Ee" 
J = S+p. 


At time T, 
I= J =max{Sr, E} + De"), 


This yields the put-call parity for all time. 


. When there is no dividend, we have shown that 


C-P<S—Ee tt, 
When there is dividend payment, we know that 
Co<C, Pp2=P 
Hence, 
Co-Pp<C-P<S-—Eet-%, 
Here, Cp and Pp refer to Amperican call and put with paying dividends. For the American 
call option, we should not exercise it early, because the dividend will cause the stock price 
to jump down, making the option less attractive. We should exercise it immediately after the 
ex-dividend date. 


. For the American put option, we consider 


I=C+D+EH, J=P4+S. 
If we exercise P at 7 < T after the dividend being paid, then S$; < FE and 


I(r) = De) + Ber), 
J(r) = E+De"-%, 


We have J(7) < I(r). Hence J(t) < I(t) for all t < 7. 


. If the put option is exercised before the dividend being paid, then we should consider [ = 


C+EandJ=P+S. Att, 

I(r) = Ee ~), 

Ta) = Be, 

Again, we have J(rT) < I(r). Hence J(t) < I(t) for allt < 7. 
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3.3. The Black-Scholes Equation 


3.3.1 Black-Scholes Equation 


The fundamental hypothesis of the Black-Scholes analysis is that there is no arbitrage opportunities. 
Besides, we make the following additional assumptions: 


(1) The asset price follows the log-normal distribution. That is, the asset price S' satisfies the 
s.d.e. 
dS = uSdt + o0Sdz. 


(2) There exists a risk-free interest rate r. 
(3) No transaction costs. 

(4) No dividend paid. 

(5) Shorting selling is permitted. 


Our purpose is to value the price of an option (call or put). Let V(S, t) denotes for the price of an 
option. The randomness of V(S(t),t) would be fully correlated to S(t). Thus, we consider a port- 
folio which contains only S and V, but in opposite position in order to cancel out the randomness. 
Then this portfolio becomes deterministic. To be more precise, let the portfolio be 


== As. 


that is, long an option and short A amount of the underlying asset. In a small time step dt, the 
change of this portfolio is 
dll = dV — AdsS. 


Here A is held fixed during the time step. From It6’s lemma 


1 
dll = Vidt+VsdS + 5Vsso°dt — Ads 


OV 


OV OV 1 OV 
— — pSA ie t. 
+ (us 5 HOOT ppg age) 
Now, we can eliminate the randomness by choosing 
OV 
A=-— 3.10 


at the starting time of each time step. The resulting portfolio 


T=VeVvs 
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has the property that its change 


_ f OV 1 kp OPV. 
am = (4 5025?) at 


is wholly deterministic. From the hypothesis of no arbitrage opportunities, the return, a should 
be the same as investing IT in a riskless bank with interest rate r, i.e. 


- = rdt. 


Otherwise, there would be either a net loss or an arbitrage opportunity. Hence we must have 


2 
rlldt = (use usd 4 cl + 1 2.429 =) dt 


Os Oe #2 Os? 
2 AO gee V 
= (F +59 S gr) tt 
or 
OV 1 anV _ av 


This is the Black-Scholes partial differential equation (P.D.E.) for option pricing. 
e Its left-hand side V; + es Vg is the return from the hedged portfolio per unit time, 
e while its right-hand side r(V — SVs) is the return from bank deposit per unit time. 
Note that the equation is independent of pu. 
Remark. Notice that the Black-Scholes equation is invariant under the change of variable S +> AS. 
3.3.2 Boundary and Final condition for European options 
e Final condition: 
c(S,T) = max{S — E,0} 
p(S,T) = max{E — S,0}. 


In general, the final condition is 
V(S,T) = ACS), 


where A is the payoff function. These final conditions are derived from the “no-arbitrage 
opportunities” assumption. 


e Boundary conditions: 


Gi) On S =0: 
C07) 0, Nese: 


This means that you wouldn’t want to buy a right whose underlying asset costs nothing. 
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(ii) On S =0: 
p(0,7) = Ee~"(F-7), 


This follows from the put-call parity and c(0,t) = 0. 
(iii) For call option, at S = oo: 
c(S,t) ~ S— Ee"? 9), as § 00. 
Since S' — ox, the call option must be exercised, and the price of the option must be 
closed to S — Ee~"(7-"), 


(iv) For put option, at S' = oo: 
p(S,t) + 0, as S 4 co 


As S' —> oo, the payoff function A = max{ FE — S,0} is zero. Thus, the put option is 
unlikely to be exercised. Hence p(S,T) — 0as S > oo. 


3.4 Exact solution for the B-S equation for European options 


3.4.1 Reduction to parabolic equation with constant coefficients 


Let us recall the Black-Scholes equation 


OV di. a3 OPV OV 
S S V=0, (0,7 0, co). 3.12 
a 2°" ag ag" ee) ae 
This P.D.E. is a parabolic equation with variable coefficients. Notice that this equation is invariant 
under S — AS. That is, it is homogeneous in S with degree 0. We therefore make the following 


change-of-variable: 


t= log 5. 


The fraction S/F makes x dimensionless. The domain S' € (0,00) becomes x € (—oo, 00) and 


av _ asev _ av 
Ox 0x2 OS OS’ 


av _ a av 
0x2 CO: OSs 


OS OV OS OV 


Hn OS” Oe OS? 
OV oo PV 
OV pV 


Next, let us reverse the time by letting 
T=T-t. 
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Then the Black-Scholes equation becomes 


OV. Dg OV ¢ 1 oe 


Or 2° Ox? 2° Ox Me 


We can also make V dimensionless by setting v = V/E. Then v satisfies 


2 
Le ¢ 50°) OU ay (3.13) 


Or 2 - Ox? 2 Ox 


The initial and boundary conditions for v becomes 
v(x, 0) = A(a) = A(Be*)/E, 


“TT as © > OO} 


e call option: A(z) = max{e* — 1,0}, v(—oo,T) = 0, v(a,T) + e* —e 
e put option: A(z) = max{1 — e”, 0}, v(—oo, rT) = e7", u(z,T) 4 Oas x > oo. 


Our goal is to solve v forO <7 < T. 


3.4.2 Further reduction 


Let us divide the v equation by 07/2: 


200 Ov (3 ) Ov 2r 


CE 


= T 
o2 OT ~=Ox? a2 Ox oc 


Let us make a change of variable and define the following new constants: 


1 9 
ea Too ) 
a=1-—r/(=07) 
1 
b=r/(50°) 
Then the equation becomes 
Us + QVz + DU = Ver. (3.14) 


The part, v, + av,z is call the advection term. The term bv is called the source term, and the term 
Ue 1S called the diffusion term. The advection part: 


Ur + ave = (0; + ad,) 


is a direction derivative along the curve (called characteristic curve) 


dx 
=a. 


i 
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This suggests the following change-of-variable: 


y = x£-—as 
s = 5s. 
Then the direction derivative become 
Os = O,t+adz 
Oy = Oy 


Hence the equation is reduced to 
Us + bu = Vyy. 


Since s’ = s, we will still use s below instead of s’. Therefore, the equation becomes 
Us + bu = Vyy. 


Next, the equation v, + bu suggests that v behaves like e?® along the characteristic curves. Thus, it 
is natural to make the following change-of-variable 


Then 
du = 0,(e%v) = eA, + bev = e°5(8, + b)u = e*O2u = Ku. 


Thus, the equation is reduced to 
Ue = thas: UE (=00,00) 4 > 0. 


This is the standard heat equation on the real line. Its solution can be expressed as 


(oe) il _z)2 
a= / oa ID, asd, 


—oo VAs 


where f is the initial data. 


3.4.3 Black-Scholes formula 
Let us return to the Black-Scholes equation (5.6). Thus, the solution v 


(,7) =e / (eee eC (3.15) 

v(“z,T) =e "7 ———e 2027 z) dz : 
00. V/ 2re*r 

In terms of the original variables, with the change of variable S’ = Ee*, we have the following 

Black-Scholes formula: 


(nS, )-+(r— $0? )(P—4))? 


202(T—t) A(S") ds’ (3.16) 


V(S,t) = 


err | 1 7 
0 /2no0?(T — t)S" 
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We may express it as 


V(S,t) =e 77-9 if P(S',T, S,t)A(S") dS’. (3.17) 
0 


Here, 


PS t, 54) = (3.18) 


1 oxy (LM B=O = 30°F - HP 
2no?(T — t)S! : 207(T — t) 


This is the transition probability density of an asset price model with growth rate r and volatility 
a. In other words, V is the present value (deducted by e~ (7-9) of the expectation of the payoff A 
under an asset price model with volatility 7 and growth rate r. It is important to note that it depends 
on 7, not on the growth rate of the underlying asset. We shall come back to this point later. 


3.4.4 Special cases 


1. European call option. The rescaled payoff function for a European call option is 


A(z) = max{e? — 1, 0}. 
Then 


OO 1 12 2 2 
ren ng (a@-2-(30%-1r)1)*/ (2077) (2 
u(xz,T) =e e 2 e 1 dz. 3.19 

( ) | V 21027 


This can be integrated. Finally, we get the exact solution for the European call option 


c(S,t) = SN(d,) — Ee" 9 N (dp), (3.20) 
where 
N oo i =24 3.21 
(y) —_ ae 2 zy ( : ) 
= log(3) + (r + $07)(T - t) 
d = ee . (3.22) 
log(§) + (r — 402)(T -t) 
a : 2 
2 oVT—t ere 


The calculation is done as below. Let us use the following notations for abbreviation. 
a=r—o'/2, D=0"7, dVD=x+ar+D, d&VD=x+ar 


We have 


u(a,r) = fo : exp ( pl? £ ar)?) (= 1)dz = 1-0. 
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= i D? + 2D(a + ar) 
= tar +D))? d 
/ =p exp ( SD (z—(a+ar yr + xD 2 


-[- po ( ame in VD) + (dh + dy)VD/2) dz 
exp ( pl? anvD)) ay 


— e(titde)VD/2 i? 
TT ae 1 2 
—— exp —5(2/VD - dh) dz/VD 
0 


1 
0 vV27D 


o 


S 
# 
S 
E 
n= [f° : exp ( 
0 27D 


= a exp ( sp ‘nvD}) ee 
= | eo (-5e/vD- a) dz/VD 
Z [ Sg ie Ne, 


—dy V2 


V(S,t) = Ee~"7-9 (I — ID) = SN(d1) — Ee" ON (dy). 


Remarks 
e We have seen that a lower bound of cis c > S — Ee7t(f-?), 


e The formula is a modification of the previous lower bound formula. The function NV is 
the error function. (+00) = 1. The two parameters d; — dz = oT — t, the standard 
deviation of the asset price model. 


Exercise. Show that 


(a) N(0) = 1/2, N(x) > las 2 > oo; 
(b) N(x) +N(—x) =1foranyx eR 


Exercise. Prove the formula (3.20). 
2. European put option. Recall the put-call parity 


c+ Ee?) = p+ 8. 
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We can obtain the price for p from c: 


p(S,t) = Ee-"?-) N(—dy) — SN (—d)). (3.24) 


Exercise. Prove (3.24). 


. Forward contract Recall that a forward contract is an agreement between two parties to buy 


or sell an asset at certain time in the future for certain price. Its value V also satisfies the B-S 
equation. The payoff function for such a forward contract is 


A(S)=S—E. 


In terms of the rescaled variable, A(z) := A(Ee*)/E = e*—1. Letz = In(S/E),7 = T-t) 
and u(xz,7) = e’’V(Ee*,t)/E. Then u satisfies the advect heat equation. Its solution with 
initial data A(z). Its solution is given by 


iL CO _ (@=a-(r=9? /2)r)? 
ula, 7) = Jone e 2027 (e = 1) dz 
= 6 24, 
Hence, 
V(S,t) =S— Ee ?F-%), (3.25) 


This means that the current value of a forward contract is nothing but the difference of S and 
the discounted . Notice that this value is independent of the volatility o of the underlying 
asset. 


Exercise. Show that the payoff function of a portfolio c — pis S — E. From this and the 
Black-Scholes formula (3.16), show the formula of the put-call parity. 


. Cash-or-nothing. A contact with cash-or-nothing is just like a bet. If Sp > E, then the 


reward is B. Otherwise, you get nothing. The payoff function is 


B ifS>E 
A(5) = { 0 otherwise. 


Using the Black-Scholes formula (3.16), we obtain the value of a cash-or-nothing contract to 
be 
V(S,t) = Be"? -9N (da). (3.26) 


Exercise. Verify this formula. 


5. Supershare. Supershare is a binary option whose payoff function is defined to be 


_[ B if, <S<E 
A(5) = { 0 otherwise. 
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One can show that the value for this binary option is 
V(S,t) = Be"? (N (do(E1)) — N (do(E2))) , 
where d2(£) is given by (3.23). 


Exercise. Verify this formula. 
6. Deterministic case (o = 0). In this case, the Black-Scholes equation is reduced to 
Vi.trSVs —rV =0. 
Or in 7, x and wu variables: 
Ur — TU, = O 


with initial data 
u(x,0) = A(Ee”)/E, 


Thus, its solution is given by 
u(x,T) = u(x + 17,0) = A(Ee*t"")/E = A(Se")/E. 


Or 
V(S,t) =e 774 A(Se"F-9), 


This means that when the process is deterministic, the value of the option is the payoff func- 


tion evaluated at the future price of S at T (that is Se"(7—")), and then discounted by the factor 
eT (T-t) : 


3.5 Risk Neutrality 


Notice that the growth rate js of the underlying asset does not appear in the Black-Scholes equation. 
The option may be valued as if all random walks involved are risk neutral. This means that the drift 
term (growth rate) js in the asset pricing model can be replaced by r. The option is then valued 
by calculating the present value of its expected return at expiry. Recall the lognormal probability 
density function with growth rate r, volatility o is 


1 [in($) — (r — 50°)(F -)P 
PS T.84) = ex 2 2 ; (3.27) 
) 2na?(T — t)S! ( 207(T —t) 
This is the transition probability density of an asset price model in a risk-neutral world: 
dS 
at rdt + odz. (3.28) 


The expected return at time 7 in this risk-neutral world is 


/ P(S',T, S,t)A(S')dS". 
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At time t, this value should be discounted by ets) 
V(S,t) =e 77-9) / P(S',T, S,t)A(S")ds". 


We may reinvestigate the function VV and the parameters d; in the Black-Scholes formula. After 
some calculation, we find 


N (do) = [ rez, S,t)ds'. (3.29) 
E 


This is the probability of the event {S > E}, where S obeys the risk-neutral pricing model: 


=rdt+adz. 


cn oe 


Similarly, one can show that 


S',T, S,t)S'as! 
Ser(T-t) 


N (di) = La Bl (3.30) 


is the expectation of S at T when S = 1 at t and under the condition that S > E at T. 


Exercise. Check formulae (3.29) and (3.30). 


3.6 Hedging 


Hedging is the reduction of sensitivity of a portfolio to the movement of the underlying of asset 
by taking opposite position in different financial instruments. The Black-Scholes analysis is a dy- 
namical hedging strategy. The delta hedge is instantaneously risk free. It requires a continuous 
rebalancing of the portfolio and the ratio of the holdings in the asset and the derivative product. The 
delta (A) for a whole portfolio is A = a This is the sensitivity of II against the change of S. By 
taking dII — A- dS = 0, the sensitivity of the portfolio to the asset price change is instantaneously 
zero. 
Besides the delta hedge, there are more sophisticated hedging strategies such as: 


e711 
Gamma: T = a5? 
Theta:O = re 
Vega: = 

II 

tho:p = as 


Hedging against any of these dependencies requires the use of another option as well as the asset 
itself. With a suitable balance of the underlying asset and other derivatives, hedgers can eliminate 
the short-term dependence of the portfolio on the movement in t, S, a, r. 
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3.6.1 The delta hedging 


Suppose a financial institution write 1,000 European calls on 10,000 shares of a stock. How does 
this institution hedge the risk from the price fluctuation of the underlying stock? The Delta-hedge 
is such a strategy. It is a dynamic hedging strategy. The institution long A shares of the underlying 


stock. The A is chosen to be 


Oc 
A= 55° 


If A = 0.5, this means that it should long 0.5 x 10,000 = 5,000 shares. In this case, the total 
portfolio II = c— AS is instantaneously neutral. That is d1[ = dc— AdS' = 0 under a small change 
of dS. We call IT is delta neutral. However, this delta neutral holds only for a short period of time. 
It should be adjusted periodically. This is called a dynamic hedging. 

For the Delta-hedge for the European call and put options, we have the following propositions. 


Proposition 1. For European call options, its A hedge is given by 
A =N(di). 


Proof. By definition, A = 0c/OS. Since c = SN(d,) — Ee~"'7-N (dz), we get 


oc = N (di) +S. N" (dy) . dig — Ee7"(T-9 N" (da) dag. 
Since : . 
1, — ee +Ct+ SVE -2) _ log(S) + (7 - ST - 4) 
: oVT —-t , oVT —t 
we get 
dig = dg : » Nd) = te GP 
So/(T —t) V2n 
Hence, 
ve = Ndi) + (SN(d,) — Be" ?-ON(da)) [(SoVT 1) 


= N(d\)+Z/(SovT — t). 
We claim that Z = 0. Or equivalently, 
S N"(d1) 


EN'(d2) 


=r 


This follows from the computation below. 


S N"(d;) 
E N'(dz) 


= ew eo (dj 43) /2_ 
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From @22), 23) 


1 
@-d = — ((2 +r74 aT) (w+rr >r)?) 
= 2(x+r1r7T) 
Hence, 
SN'(di) Ie a ee | 
E N'(d2) 
Proposition 2. For European put options, its A hedge is given by 
A = N(-d1). 
Proof. From the put-call parity, 
Op Oc 
A= == —5-1=N(d,}) -1=—N(-d 
as as N( 1) N( 1) 


3.7 Time-Dependent r, o and for Black-Sholes equation 


We can extend Black-Schole analysis to the case when r, o, jz are functions of t, but still determin- 
istic. We use the change-of-variables: 


S=Ee*, V=Ev, tT=T-t. 


The Black-Scholes equation is converted to 


2 2 
= sd (7) Vex + (r(T) e {)),, r(r)v (3.31) 
We look for a new time variable 7 such that 
d? = 0*(r)dr 
For instance, we can choose 
[= | o*(r) dr 
0 
Then the equation becomes 
1 
Us = 9 Vex + a(7)vz = b(7)v. (3.32) 


To handle the term a(7), we consider the following characteristic equation: 


dx : 
a —a(T) 
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This can be integrated and yields 


r= -{ a(r’)dr’ + y, 
0 


where y is an integration constant. The lines with constant y’s are called characteristic curves. 
Along the characteristic curves, the derivative 0; + a(7)0, is just the derivative in time with fixed 
y. In fact, we introduce the following change-of-variable: 


=] cat H with { yout na a(r!)dr! = 2+ A(#) 


T T, =T 
Then, 
0 fees Oy OO 
Ox" Ox Oy Oy’ 
and 
Q, _ Fa Ox O O 0 


Bal gn Or Ones ro ae 
The equation (3.32)is transformed to 
i 
Ux = =Uyy — O71 )y. 
2 
Let B(71) = 7" 


9. 0(7')dr’ and u = e? (7)y, then us, = 4uy,. And we can solve this heat equation 
explicitly. 


2 


3.8 Trading strategy involving options 


An investor can design his/her payoff by using options and underlying stocks. The options whose 
payoff are max{ Sp — £0} or max{E — Sr, 0} are called vanilla option. They are the simplest 
payoff functions. In this section, we shall discuss how to design more general payoff functions to 
fit an investor’s anticipation. 


3.8.1 Strategies involving a single option and stock 
There are four ways: 


a. II = S — c (writing a covered call option). When an investor short a call, we say he write 
a nake call. The investor has cash inflow at the beginning. However, short selling is highly 
risky. If the stock rises S; > E, then the investor needs to buy the share on S'r and sell to the 
call holder on E.. The investor loses S7 — EF’. This loss can be unlimited because S7 may rise 
very high. To limit this risk, the investor can long a share S; at time ¢. That is, his portfolio 
is I = —c+ S. At expiry, if Sp > E, then he can use this share to sell to the call owner on 
F, instead of Sp. His risk is limited. We say that this share covers the call, and say the writer 
writes a covered call option. The payoff of Lis A = Sp — max(Sr — E,0) = min{ Sr, E}. 
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b. II = c— S (reverse of a covered call). Suppose an investor anticipates the stock price will 
decrease. So he shorts a share and receive money 5; at time t. If S77 does decrease, he earns 
money. If unfortunately, the stock price increases, then he loses Sr — S; at time T’. This risk 
is unlimited. Thus, short selling has unlimited risk. However, he can buy a call to limit this 
risk. That is, his portfolio II = —S +c. Since c < S, he still receive some money from II at 
beginning. At expiry, if the price goes beyond EF, then he can exercise this call to cover the 
shorted share. The payoff now is A = — min{.S, E} > —E. Thus, his loss is at most EF. 


c. Il = p+-S (protective put). When an investor anticipates a stock will increase, he set up a 
portfolio II = S. This costs him S;. If the prices Sp goes up, he earns Sy — S;. However, 
if the price goes down, he sells it and loses S; — Sr. If he wants to limit his loss, he can set 
up a portfolio II = S + p. This costs him more at beginning. But at expiry, if the price goes 
down, he can sell the share on F instead of a lower price Sy. Thus, he limits his loss. In this 
portfolio, the payoff is A = S + max{E — S,0} = max{S, E}. The downside is limited (by 
£) and the upside is unlimited (Sp can be arbitrary large). 


d. Il = —p — S (reverse of a protective put). In this portfolio, the investor anticipates the stock 
price will decrease. He shorts a share S; and also a put. He receives money from shorting 
these. At expiry, the payoff is — max{S, E}. This means that he pays E if Sp < E. However, 
if S'p is large, then his risk is unlimited. 


Below are the payoff functions for the above four cases. 


3.8.2 Bull spreads 


In this strategy, an investor anticipates the stock price will increase. However, he would like to give 
up some of his right if the price goes beyond certain price, say Fz. Indeed, he does not anticipate 
the stock price will increase beyond EF. Therefore he does not want to own a right beyond Ez. 
Such a portfolio can be designed as 


Il = Cp, — Cr, Ey < Eo, 


where Cp, is a European call option with exercise price E; and C'z,, Cz, have the same expiry. 
The payoff 


A = max{Srp — F,0} — max{Sp — E>, 0} 


0 if Sp < Ey 
S;p — FE if Ey < Sp < E> 
Fo—-E, if Sp > E> 


Since FE, < E2, we have Cp, > Cz. A bull spread, when created from Cz, — C'g,, requires an 
initial investment. We can describe the strategy by saying that the investor has a call option with a 
strike price E, and has chosen to give up some upside potential by selling a call option with strike 
price Ey > EF. In return, the investor gets Hy — EF if the price goes up beyond Fo. 
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(a) (b) 


(c) (d) 


Example: Cg, = 3, Cz, = 1 and £) = 30, Eg = 35. The cost of the strategy is 2. The payoff is 


0 if Sp < 30 
Sp —30 if 80 < Sp < 35 
5 if Sp > 35 


The bull spread can also be created by using put options 


l= Pr, — Pr, Ey < Es. 


3.8.3 Bear spreads 
An investor entering into a bull spread is hoping that the stock price will increase. By contrast, an 
investor entering into a bear spread is expecting the stock price will go down. The bear spread is 


Il = Cr, _ Cr, Ey < Eo. 


There is cash flow entered (C'z, — Cp, ). The payoff is 
A =—max{Srp — £;,0} + max{Srp — E>, 0} 
0 if Sp < Ey 
— —-SrptF, iff, < Sr < ko 
—Fo+ EF, if Sr > Eo 
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3.8.4 Butterfly spread 


If an investor anticipate the stock price will stay in certain region, say, Ey < Sp < Es, he or she 
can have a butterfly spread such that the payoff function is positive in that region and he or she gives 
up the return outside that region. 


1. Butterfly spread using calls: Define the portfolio: 
= Cr, = 2C'Ry + Crs, with Ey < Eo < Ez. 


where £3 = Ey + (E2 — E)). Its payoff function is a piecewise linear function and is 
determined by A(E,) = A(£3) = 0, A( £2) = E2 — E;. Below is the graph of its payoff 
function. 
E,-E, 
S-E, E,-S 
E E, Ey Ss 


Example: Suppose a certain stock is currently worth 61. A investor who feels that it is 
unlikely that there will be significant price move in the next 6 month. Suppose the market of 
6 month calls are 


Lace 
55 10 
60 7 


65 5 
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55 60 65 Ss 


The investor creates a butterfly spread by 
Il = Cr, — 2C rg, + Cr,. 
The cost is 10 + 5 — 2 x 7 = 1. The payoff is the figure below (Figure (1p). 


2. Butterfly spread using puts. 


Ey, + E: 
Pr, + Pe,—2Pp,, FE, < £3, Eg= Sa 
linear 
0 if S < By —AE, or S > By —-AE 


Remark 1. Suppose European options were available for every possible strike price F, then any 
payoff function could be created theoretically: 


A(S) = S> Aiw(S — Ej) 
where £; = iAE, A; is constant and 


1 if S— BE; =0 
y(S-E;)=% 0 if |S — Bj) > AE 
linear for |S — E,| < AE. 


Then A(F;) = A; and A is linear on every interval (£;, F;,1) and A is continuous. As AE — 0, 
we can approximate any payoff function by using butterfly spreads. 
Remark 2. One can also use cash-or-nothing to create any payoff function: 


A(S) = 5° Avh(S — Ei), 
where 


1 if0<S<AE 
0 otherwise. 


(8) = H(S) — H(S — AE) = { 
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and ff is called the Heaviside function. The value for such a portfolio is 


Vo = ent | p(8',7, S,2)A(S)dS’ 
eT TODA P(E; < S < Ej41). 


3.9 Derivation of heat equation and its exact solution 


3.9.1 Derivation of the heat equation on R 


The heat equation on R reads 


{ Ut =Ure, « ER,t>O, 
u(x, 0) = f(x). 


This equation models heat propagation on the line. Here, u(-, t) represents the temperature distribu- 
tion at time ¢ and f the initial temperature distribution. The derivation of this equation is based on 
a physical law — the conservation of energy, and the Fourier law of heat flux. Given a temperature 
distribution u(x,t), the Fourier law describes that the temperature flows from high temperature to 
low temperature at a rate q called heat flux. It is the energy passing through a point (per unit area) 
per unit time. The Fourier law states that 


q = —KUzy. 
Here, the constant « is called thermal conductivity. It is a positive constant. It differs for different 
material. For instant, « of cooper is much larger than that of woods. According to thermal dynamics, 
the energy density e is related to temperature by e = c,u. The constant c, is called specific heat (at 


constant volume). The conservation of energy states that the change of energy in an interval (a, b) 
per unit time is the same as the heat flux flows into (a, b) from boundaries. Mathematically, it is 


a [ ear = ala ae - [gece 


b b 
| cvdmae = | (KUg)e dx. 


This holds for any interval (a, b). Thus, the integrands much be the same. This gives 


This gives 


Ou = Dury. 


where D = k/cy > 0 is the heat conductivity. 
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3.9.2 Exact solution of heat equation on R. 


If we rescale time by t’ = Dt, then the heat equation becomes 
OU = Ue: 


Thus, without loss of generality, we consider the heat equation with D = 1. We shall show that its 
solution is given by 


u(at)= f Ge-wtfu)dy, Gle,t) = Foe ws, 


The function G(-,t) is the Gaussian distribution. It corresponds to the temperature distribution 
at time t when the initial temperature distribution is a hot spot at x = 0. The derivation of this 
formulation is decomposed into the following steps. 


1. Superposition principle: If u and v are solutions, so is au + bv, where a and b are constants. 
This is called superposition principle for linear partial differential equations. 
Now, we imagine f can be approximated by piecewise step function: let h be a small mesh 
size, 7; = jh and f can be approximated by f), defined by 


fein, frlo 2 ane eh 


j=-oo 


Here, 


1 1 for —h/2<a<h/2 
Gi) Oy SEAN Gems 


Suppose each solution corresponding to f(x;)dp(a) is f(xj;)Gp(a — xj). Because the equa- 
tion is linear, then the general solution is just the linear superposition of these solutions: 


oo he f G(a — y, t) f(y) dy. 


Here, G),(“—«;,t) is the solution of the heat equation corresponding to the initial data 5), (a — 
x;). Notice that we have used another important of heat equation, namely, the translation 
invariance: if u(x, t) is a solution of the heat equation with initial data u(x, 0), then u(x—a, t) 
is a solution of heat equation with initial data u(a — a,0). This is because the equation is 
unchanged under the translation transformation x + x — a. Now if G;,(z,t) is the solution 
with initial data 5;,(x), then Gp,(a — 2;,t) is the solution corresponding to the initial data 
Of Dap): 


2. Translational invariance. Because the equation is translationally invariant, we can solve the 
equation with initial data to be y;,(z)/h. This initial function tends so-called 5-function as 


52 


CHAPTER 3. BLACK-SCHOLES ANALYSIS 


h — 0. Thus, what we need is to solve the case when the initial function is a 6 function. A 
6-function is a generalized function having the property: 


ae) = { 0 «#0 


co r=0 


[ @ae= 


The 6-function is considered to be the limit of 5;,(a) as h — 0+ in the sense that 


[_ He)(x) ae = fim J Fla) bn(x) de = 5 raz fi fed = $(0), 


for any smooth function f. 


. Parabolic scaling. To solve the heat equation with an initial hot spot at x = 0, we use 


the parabolic invariance property of the heat equation. Namely, the equation and the initial 
condition are both invariant under the transformation 


Bi Ney EES NE. XO: 


This suggests that the solution is a function of € = a/,/t. So, let us look a solution of the 
form u(x,t) = U(a/,/t). Plug it into the equation, we obtain 


This leads to 


-2u'(9 =U") 
(InU’(€))' = = 
Integrate it, 
2 
InU’(€) = < +Const., U'(é)= Car 


Hence Z 
= if Cer /4 dn. 
—Co 


Such a solution has the property: U(—oo) = 0 and U(+o0) = CvW4. We choose C = 
1// 4x for normalization. However, this is still not what we want. But we notice that if 
v(x, t) = U(a/Vt) is a solution, so is v,. Thus, we let 


G(a,t) = OnU (=) = = aa x? /(At) 


Then G is a solution and satisfies 
CO 
/ Ge de= 
—Co 


3.9. DERIVATION OF HEAT EQUATION AND ITS EXACT SOLUTION 53 


4. Verification. Finally, one check 


= [™, G(x — y,t) f(y) dy is a solution of the heat equation 


e u(x,t) > f(x) as t > 0+ for any bounded continuous function f. 


To show the first one, we check that G; = G',,. This is left to you. 
For the second one, the integral 


ie Bett fy)ay= fot G(-y, t+ ydy =f Gly,t)fle+y) dy 


represents the average of f about x with weight Gy, t). We notice that G(y, t) is a Gaussian 
distribution with standard deviation /t. As t > 0, this weight is more and more concentrated 
at 0. Therefore, this average becomes the average of f roughly in a \/t-neighborhood of z. 
Its limit is f(a) for any bounded smooth function /f. 


We can also give a more rigorous proof. We want to show 


lim Ju(e,*) — rei=|[~ G(x —y,t)fly) dy — f(@)| = 


Notice that = 
/ G(a — y,t) dy = 1foralla € R,t > 0. 
—oo 
Hence 


[ ee@-vtway- te) = [ ee-y UW Fe) aw. 


—co 


Therefore, os 
lu(e,t) — f(e)| < / G(e — 9,8) |f(y) — F(@)| ay. 


Here, we have used G(x, t) > 0. To estimate this integral, we break it into two parts 


io G(x—y,t) |f(y) — f(x)| dy = (/ , Ea | f G(a—y, t) |f(y) — f(«)| dy =I+I, 
=e L—Yl< r—y|> 


where 6 > 0 is a small parameter to be chosen later. Since f is continuous, we have for any 
€ > 0, there exists a 6 > 0 such that 


lf(y) — f(x)| < €, whenever |x — y| < 6. 


Thus, 


r= | Gle— w IFW) - fla)ldy <e | Gla —y,t) dy <. 
|ja—y|<d |ja—y|<d 
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For II, since we assume f is bounded, there exists a constant MM > 0 such that | f(a)—f(y)| < 
2M. Thus, 


T<2M G(a — y,t)dy =2M G(y, t) dy. 
|x—y|>6 ly|>6 


This integral is the tail estimate of a Gaussian distribution: 


1 2 2 
G(y,t) dy = 2 | e VAD dy = —=e* dz. 
(yt) y>o VArt z>6/(2Vt) Va 


As t + 0+, this integral goes to 0. Therefore, we can choose a small t; > 0 such that 


ly|>d 


I <eforallO <t< fy. 


Combining J and IT, we have for any € > 0, there exists £; > 0 such that for any 0 < t < ty 
we have J + II < 2c. This shows J + IT > Oast > 0. 


Chapter 4 


Variations on Black-Scholes models 


4.1 Options on dividend-paying assets 


Dividends are payments to the shareholders out of the profits made by the company. We will con- 
sider two “deterministic” models for dividend. One has constant dividend yield. The other has 
discrete dividend payments. 


4.1.1 Constant dividend yield 


Suppose that in a short time dt, the underlying asset pays out a dividend DoSdt, where Do is a 
constant, called the dividend yield. This continuous dividend structure is a good model for index 
options and for short-dated currency options. In the latter case, Do = ry, the foreign interest rate. 

As the dividend is paid, the return a must fall by the amount of the dividend payment Dodt. It 
follows the s.d.e. for the asset price is 


d 
< = ( — Do)dt + odz. 


To value the corresponding option V, we consider the portfolio : II = V — AS as before. In one 
time step, the change of the portfolio is 


dil = dV — AdS — AD) Sdt, 
where the last term —A Dp Sdt is the dividend our assets received|'| Thus 


dll 


dV — A(dS + DoSdt) 
os (55 _ A) dz 


OoVeoi OV 
@ Do) 57a + 57 5 aS? + Vi (ps Do)AS = ADoS) dt 


Te you own a share, after dt, you will receive S Dodt dividend. 
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1 2: 
= (vi ae 9222 at 


as’? 22 OS? 
Here, we have chosen A = a to eliminate the random term. From the absence of arbitrage 
opportunities, we must have 
dil = rIldt. 
Thus, 
OV 1 OPV OV 
V, — DoSaa + =0°S° =r (V—-Se)-. 
Be BI. Pee 5 G2, ( = | 
1.e., 
1 OPV ov 
V+ =07S? | Do)S V=0 
bg © gga PP ag” 


This is the Black-Scholes equation when there is a continuous dividend payment. 
The boundary conditions are: 


e(O;t); = 0, 
c(S,t) ~ Se Pott) 
The latter is the asset price S discounted by e~?0(7—*) 


payoff function c(S,T) = A(S) = max{S — E, 0}. 
To find the solution, let us consider 


from the payment of the dividend. The 


c(S,t) = e POF 9 6, (S, 2). 


Then c, satisfies the original Black-Scholes equation with r replaced by r — Do and the same final 
condition. The boundary conditions for cy are 


ca(S,t) ~ SasS— oo 


Hence, 
(S,t) = SN(d19) — Ee~°-P0) 9 N (dy 0) 
where 
oe In 2 + (r — Do + 507)(T — t) 
; oVT —-t : 
dz = dio —ovT —-t. 
The original 


c(S, t) = Se~PolF-D N (dy 9) a Ee?) N(dz9). 


This means that the call option on a dividend-paying asset is just like a call option in risk-free 
environment with rate r — Do and discounted by e~"(7—*), 
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Remark. We notice that c \, as Do 7. This is natural because paying dividend will lower the 
value of S', hence c. If the dividend is paid more, the value of call becomes lower. 


Exercise. Derive the put-call parity for the European options on dividend-paying assets. 


4.1.2 Discrete dividend payments 


Suppose our asset pays just one dividend during the life time of the option, say at time tg. At tg+, 
the asset holder receiver a payment d,S(tg—). Hence, 


S(tat+) = S(ta—) — dyS(ta—) = (1 — dy) S(ta-). 
We claim that across the jumps, V should be continuous, i.e., 
V(S(ta—), ta—) = V(S(tat), tat). 


Otherwise, there is a net loss or gain from buying V before ¢, then sell it right after tg. To find 
V(S, t), here is a procedure. 


1. Solve the Black-Scholes from T to tg+ to obtain V(S, ta+) (using the payoff function A) 


2. Adjusting V by 
V(S, ta—) — Lace ry dy) 8S, tat) 


3. Solve Black-Scholes equation from tq to t with the final condition V ((1 — dy)S,tg+). 


For the case of call option, we can find its explicit expression. Let cg be the European option for 
this dividend-paying asset. Then 


ca(S,t) = c(S,t, E) fortg+ <t<T 
calS, ta=) = eg(S0— dy) tat) =e( SU = d,) ta, 2) 
Note that 
e(S(1 — dy), T, E) = max{S(1— dy) — E,0} = (1 — dy) max{S — (1— dy)" E, 0} 
and the linearity of the Black-Scholes equation, we obtain 
c(S(1 — dy), tg, Z) = (1 — dy)c(S, tg, (1 — dy)" E). 


Therefore, 
ca(S,t) = (1 — dy)e(S,t, (1 — dy)“ B). 


That is, the call option price is reduced to its original value with strike price E’ replaced by E'/(1 — 
dy). 
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4.2 Futures and futures options 


4.2.1 Forward contracts 


Recall that a forward contract is an agreement between two parties to buy or sell an underlying asset 
on a certain price F at a certain future time 7’. Here, F is called the delivery price. The payoff 
function for this forward contract is A = S; — E. Based on the no arbitrage opportunity, we have 
seen that the price for this forward contract is 


f=S-Eetr, 


Definition 4.3. The forward price F for a forward contract is defined to be the delivery price E 
which would make that contract have zero value, i.e., 


fr=0, or F=Se"F-, 


One can take another point of view. Consider a party who shorts the contract. He can borrow an 
amount of money 5; at time ¢ to buy an asset and use it to close his short position at T. The money 
he received at expiry, F’, is used to pay the loan. If no arbitrage opportunities, then 


F= Ser : 


If the delivery price F of a forward contract is not equal to the forward price F’, then the value of 


this forward contract is f = S;— Ee~"(—. This is the result we obtained in the last chapter. Thus, 


in order to have f = 0, we should take the delivery price to be S,e"(7— 


4.2.2 Futures 


Futures are very similar to the forward contracts, except they are traded in an exchange, thus, they 
are required to be standardized. This includes size, quality, price, expiry, ... etc. Let us explain the 
characters of a future by the following example. 


1. Trading future contracts 
e Suppose you call your broker to buy one July corn futures contract (5,000 bushels) on 
the Chicago Board of Trade (CBOT) at current market price. 
e The broker send this signal to traders on the floor of the exchange. 


e The trader signal this to ask other traders to sell, if no one wants to sell, the trader who 
represents you will raise the price and eventually find someone to sell 


e Confirmation: Price obtained are sent back to you. 
2. Specification of the futures: In the above example, the specification of this future is 


e Asset : quality 
e Contract size: 5,000 bushel 
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Delivery arrangement: delivery month is on December 


price quotes 


Daily price movement limits: these are specified by the exchange. 


Position limits: the maximum number of contracts that a speculator may hold. 
3. Operation of margins 


e Marking to market: Suppose an investor who contacts his or her broker on June 1, 2016, 
to buy two December 2016 gold futures contracts on New York Commodity Exchange. 
We suppose that the current future price is $400 per ounce. The contract size is $100 
ounces, the investor want to buy $200 ounces at this price. The broker will require 
the investor to deposit funds in a “margin account”. The initial margin, say is $2,000 
per contract. As the futures prices move everyday, the amount of money in the margin 
account also changes. Suppose, for example, by the end of June 1, the futures price has 
dropped from $400 to $397. The investor has a loss of $200x3=600. This balance in 
the margin account would therefore be reduced by $600. Maintaining margin needs to 
deposit. Certain account of money to keep that futures contract. 


e Maintenance margin: To insure the balance in the margin account never becomes nega- 
tive, a maintenance margin, which is usually lower than the initial margin, is set. 


We tabulate the futures prices, the main/loss in each day. Suppose there are n days to the 
maturity date. Let 6 = r(T — t)/n be the riskless interest rate per date. Let F; be the future 
price at the end of day 7. Suppose the margin has money / initially. The gain/loss in each 
day is the following table. 


Day 0 1 2 ine n 
futures price Fo Fy Fo see Fy-1 Fy, = Sr 
Margin (compounded) MM Me Me”? oo tee Me™ 
gain/loss 0 (Fi — Fo) (Fo — Fi) vee vee (Fn — Fn-1) 
gain/loss compounded. 

to day n Oo (RS Her ea eC os oho ie e™ 


At the end of day 2, the margin account should have 


a 
Me + $0 (Fy — Fy_sjeO*. 
k=1 


This account is required to be kept an amount money greater than certain maintenance fee. If 
it is lower than this maintenance fee, it will be closed. 


Here, the futures prices F; are the market future prices. What are the futures prices if there is no 
arbitrage opportunities? We show below that under the no-arbitrage assumption, the future price at 
time t is indeed identical to the forward price. That is F, = Sper(F-4) | 
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Theorem 4.4. Forward price and futures price are equal when the interest rates are constant. 


Proof. Suppose a futures contract lasts for n days. Let the future prices are 
Fo, eg Fy 


at the end of each business day. Let 6 be the risk-free interest rate per day. That is, d = r(T —t)/n. 
Consider the following two strategies: 


e Strategy 1: Suppose the forward price is Go. 


— We invest Go in a risk-free bond; 


— take a long position of amount e”® forward contract. 
At day n, Goe” is used to buy the underlying asset at price Spe”. The payoff is Spe”. 
e Strategy 2: Suppose the futures price at day 0 is Fo. 


— We invest fy amount of money in a risk-free bond; 


— take a long position of e*® amount of future at the end of day i — 1,i = 1,...,n. 


Day 0 1 2 (att n 
futures price Fo Fy Fy see Fy Fy 
position e e?5 e° tee er? 0 
gain/loss 0 e°(F, — Fo) e?)(Fy — Fi) vee vee ce” (F, — Fy-1) 
compound 0 &(Fr—Foje™-Y8 67 (Fy — Fi)e-7)F i... vee ec”? (Fy — Fn—1) 


The total gain/loss from the long position of the futures is 
" . 2% 
SS 0(Fi - Fi-rje®- el)? = (F, — Foe” = (Sp — Fe”. 
i=l 
At time T, we receive Fye” from the initial investment Fy. Thus, the total payoff of strategy 2 is 


Foe™ + (Sp — Fy)e™ = Spe™. 


Since both strategies have the same payoff, we conclude their initial investments must be the same, 
ie., Fy = Go = Spe"), 


4.2.3 Futures options 


Options on futures are traded in many different exchanges. They require the delivery of an underly- 
ing futures contract when exercised. When a call futures option is exercised, the holder acquires a 
long position in the underlying futures contract plus a cash amount equal to the current futures price 
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minus the exercise price. Suppose the expiry date of the call is JT}, which is T, < T’, the expiration 
of the underlying future. The payoff of the call option on future is thus 


A = max(F7, — E,0). 


In addition of this payoff, the holder also receive a future contract, yet it costs zero to enter the long 
position of a future. 


Example An investor who has a September futures call option on 25,000 pounds of copper 
with exercise price EF = 70 cents/pound. Suppose the current future price of copper for deliv- 
ery in September is 80 cents/pound. If the option is exercised, the investor received 10 cents 
x 25, 000+long position in futures contract to buy 25,000 pound of copper in September at price 
80 cents/pound. 

The maturity date of futures option is generally on, or a few days before, the earlist delivery 
date of the underlying futures contract. 

Futures options are more attractive to investors than options on the underlying assets when it is 
cheaper or more convenient to deliver futures contracts rather than the asset itself. Futures options 
are usually more liquid and involved lower transaction costs. 


4.2.4 Black-Scholes analysis on futures options 


Future Price Model As we have seen that the futures price is identical to the forward price when 


the interest rate is a constant, i.e., F = S' e’(T—-t) | From It6’s lemma, we obtain a pricing model for 
F: 


dF = (ar aos 28?) ar + eas 


0 2 05" 
= (—re™ 9) 9)dt + Se oe 
= (-rF)dt + F(pdt + odz) 
= ((u—r)F)dt + Fodz. 


Hence, 


= (u—r)dt + odz. (4.1) 


We conclude this discussion by the following proposition. 
Proposition 3. The futures price is the same as a stock paying a dividend yield at rate r. 
Next, we study the value V of a futures option. For a call option on a future, it is the right to 


buy the future on strike price EF’. Its value V is a function of fF’, t and E. We apply Black-Sholes 
analysis to value V. 
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Futures Option Price Consider a portfolio 
TI=V—AF. 
As discussed before, we choose A = a to eliminate randomness of dII. Then 


dil = dV —AdF 


OV OV 10V OV OV 
= (spHrF4 50° F*)dt + oF dz (upF'dt + oF dz) 
7 ov . 100 spay, OF OF 
~ ot” 2aF2° 


Here, jup := s — 7, is the growth rate of the undelying future. Since it costs nothing to enter into a 
future contract, the cost of setting up the above portfolio is just V.|“| Thus, based on the no arbitrage 
opportunity, 


dil = rVdt, 
and we obtain ‘, 


The payoff function for a call option is A = max{F' — E, 0}. 
To solve this equation, we recall the option price equation for stock paying dividend is 


1 2 
gg ae pee rv =0 


Ve aS 


ee OS? 


In our case, Dp = r, so the futures call option 


c(F,t) =e"? (F, t) 
where c, satisfies Black-Scholes equation with r replaced by r — r = 0. This gives 


a(F,t) = FN(d) — EN (dg), 


In — $07(T -t) 


d = ’ 
: oVI —t 
diy = dgt+tovT-t 


Notice that this option price V (i.e. c in the present case) is the same as V(.S(F,t), t), where V is 
the solution of the option corresponding to the underlying asset S. That is 


V(S,t) = SN(d\) — Ee"? -9N (da) 


es In? + (r —$07)(T -t) 
° 7 7 oVT —-t 
dj = dgt+tovT-t 


?When the futures F, = e"(7—") $4, it is the same as the forward contract price. The price f of corresponding contract 
is 0. 
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We can write this V in terms of F by S=F g Hl) Plug this into the above equation to obtain 


— In ee + (r — $07)(T —t) 


Similarly, dy = d,. Thus 


V(S(F,t),t) = SN(d,) — Ee"? ON (do) 
= Fe?) NV (d,) — Ee"? N (dz) 
= V(F,t) 


We conclude the above discussion by the following proposition. 


Proposition 4. The price for future options is the same as the price for options on the underlying 
assets. 


Finally, let us find the put-call parity for futures options. 
Proposition 5. Let c and p be the call and put options on a future. Then 
e+ Bett) =p Petr, (4.2) 
Proof. Consider two portfolios: 


A = €pRett 
B = p+Fe"?- + a futures contract 


At time JT’, 


Aa = max{Fr — £,0}+ E = max{Fr, E}, 
Ag = max{f£—- Fy,0}+F 4 (Fr —- F) = max{L, Fr}. 


Hence we obtain A = B. 
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Chapter 5 


Numerical Methods 


To evaluate the option price numerically, we need to discretize our option price model. We recall 
that there are two equivalent versions of option price models: 


e The Black-Scholes equation with initial and boundary conditions: 
Vit 28?Vas = r(V — 8V5), SS 005, 
ViGS, Fy = AGS). 
e The exact expression in terms of the transition probability of the risk-free asset price model: 
V(S,0) =e"? i; . P(S',T, S,0)A(S’)ds", (5.1) 


where P is the transition probability for the risk-free asset price model 
dS = rSdt + oSdz, 0<t<T 
S@).=S. 

That is, P(.S’, T, S, 0) is the probability density of S(T) with $(0) = S. 


The first one is a P.D.E. model. The second one is the exact solution of the P.D.E. It can also be 
viewed as a solution of a stochastic model in a risk-free world. The finite difference method is a 
numerical procedure for solving the Black-Scholes P.D.E., whereas the Monte-Carlo method is a 
numerical procedure for the second expression. The binomial method can be derived from both 
versions. 


5.1 Monte Carlo method 


We recall that the value of a European option is given by 


ViS0\=e7" / P(8,T, S,0)A(5)48 (5.2) 
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where A is the payoff function, P is the transition probability density function of S which satisfies 


= =rdt+odz, (initial state $(0) = S) (5.3) 


1.e., it is the asset price model in the risk-neutral world. The Monte Carlo simulation is a numerical 
procedure to evaluate V numerically based on the random simulation of the distribution of S(T’). 
There are two ways to find an approximation of probability distribution of S(T). 


e Use the exact solution formula for S. Recall (2.10) 
~ o2 
S(T) = Sexp (( — a F + o2(T)) ; (5.4) 


where z(7’) is the Brownian motion. We have seen that 
2(T)=VTe, €~N(0,1). 


We can simulate the distribution of S (T’) by using this formula. We generate M random 
numbers €; with distribution V’(0, 1). We approximate 


8,(T) = Sexp ((r—FyP +oVTa). 


e Solving the s.d.e. numerically. We simulate M paths (for example, M = 10,000) from 
(5.3) to obtain $;(T),i = 1,...,M. To sample M paths from 63}, we divide the interval 
(0, Z| into N subintervals with equal length At = z. We sample M7 N random numbers é., 
i=1,..,M,k =1,...,N with distribution (0,1) We then discretize by using the 
forward Euler method: 

S;(kAt) — S;((k — 1)At) 


= =rAt+odVAt, k=1,..,N,i=1,...,M 
Si((k — 1)At) . 


S;(0) = S. 
This approach certainly costs more. But it is for general s.d.e. 


The option price at time 0 is the expectation of the payoff according the probability distribution of 
S;(T), then discounted by e~"”’. Thus, 
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Remarks. 


i 


1. The error from the numerical integration by Monte-Carlo method is O ( Tu )* This can be 


proved by the central limit theorem. The error from the numerical integration of the s.d.e. 
for a sample path by using forward Euler method is O(x#): This can be proved by standard 
error analysis of numerical ordinary differential equation. Thus, the total error for Monte- 


Carlo method by using forward Euler integrator for s.d.e. is O (4 + i). If there is only 


one underlying asset, the Monte Carlo does not have any advantage. However, if there are 
many underlying assets, say more than three, the corresponding Black Scholes equation is a 
diffusion equation in high dimensions. In this case, finite difference method is very difficult 
and the Monte Carlo method wins. 


. The calibration of the volatility o of an asset from market data has been demonstrated in 
Section 


Below we provide a C code for Monte-Carlo method by using the exact formula of S (T). 


Monte_Carlo_European_option(V,S,E,r,sigma,T,M) 


{ 


} 


double St, sd, R; 

double payoff(), randn(); //randn() is the random number generator obeying normal 
distribution 

int i; 

R = (r-sigma*2 /2)*T; 

sd = sigma * sqrt(T); 


£ 
i=l; i <= M; itt) 


St = Sxexp(R+sd«randn(0,1)); 
V += payoff (St,E); 


} 
V *= exp (-r*T) /M; 


double payoff (S,E) 


{ 


} 


return max(S-E,0); 


5.2. Binomial Methods 


The binomial model is to approximate the risk-neutral asset price model dS i. S =rdt+odzbya 
binomial model, then find the option price by taking expectation of this binomial model. We first 
recall the binomial asset model. Then we illustrate how to find option price from this binomial 


model. 


5.2.1 Binomial method for asset price model 


We consider the underlying asset is risk-neutral, i,e., 


= = rdt + odz (5.5) 
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We shall approximate this continuous model by the following binomial model. 


e We choose an N, define time step size At = (T — t)/N. Define 


Ax =ovAt, 


u=e", d=e™. 


e Define a random walk: {S°, ...9} by 


n=0,..N—-1 


S™*!  ( u_ with probability p 
oS” d with probability 1 — p 


with initial price 
SCs. 


where the probability p is determined by 


erAt _d erAt = e At 


OT eed eke eae 
The reason how u, d and p are chosen in this form will be illustrated later. 
This process can be illustrated by a binomial tree. We define 
= Se" NAPS N, 


These are possible prices in every time step up to N. We start from S° = Sp = S. The price goes 
up to S; = Su = Se4*, or down to S_; = Sd = Se~*. In the second step, the possible prices 
are 

Su, Sud, Sdu, Sd?, 


which are 
S_a, So, So. 


After N steps, the possible values of the asset prices are 

{55 | |3| < N andj + N = even}. 
Each S” is a random variable. Let P;’ be the probability of the event that S” = S;: 
PP= Ps" =S)),. . |j|<1,9 + a Seven: 


Because the initial position is So, we have Py = 1. After n step random walks, the probability 
distribution of S” = S; is the binomial distribution: 


nee eS. (7)p'(1— p)”* n+j =2é 
BU by dg. = { 0 otherwise. 
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Here, ¢ is the number of steps that the price goes up in n steps. 

This binomial model depends on two parameters: u and p. ( The down ratio d = 1/u.) They 
are determined by the conditions so that the discrete model and the continuous model have the same 
mean and variance in one time step At. We recall that these conditions are 


pu+(1—p)d emt 


pu a. Gd 7 p)d? nes elarto7)At 


From these two equations, we get respectively 


From these two equations, we get 


e(2rto?)At _ G2 


erAt —d 


1 


Since d = u—“, we get a quadratic equation for u: 


1 (2r+o7)At 
ue—au+t1l=0, where a= mu 


erAt 
Solving this equation, we get 


u=5 [at vor—4]. 


From this, we recall that d and p can be expressed in terms of wu: 


1 erAt _ qd 
d=— —————— oe 
ee u—d 


Thus, p, u and d can be expressed in terms of r,o and At. 
By Taylor expansion, a simple calculation gives 


u=1+oVAt+ soot + O((At)3/2) = e7V4t 


Thus, up to O((At)~*/?) error, we can choose 


This is why we choose 
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at the beginning. The probability p is thus given by 


er At _d erAt -_ e At 
Pp ud. eAv—e-he 
A simple Taylor expansion gives 
1 1rAt 
p= en a ae 7A + O(At+ (An)? ) 
2 
We may also choose 
Big BRAG» 0 i oP’, At 
=z Z=7-+ (rT : 
D9! Reh 2 2Ax 


We should require At so that 0 < p < 1. This leads to require 


At 
— <l. 
2| Ar — 


5.2.2 Binomial method for option 


We recall that the option price is the expectation of the payoff A(S) under the risk neutral probability 
distribution P(S(T),T, So, t), then deducted by e~’(7—. This transition probability distribution 
can be approximated by a binomial distribution P;” as shown in the last section. In the discrete 
model, the risk neutral asset price only takes discrete values 5;, we shall approximate V(S;,¢ + 


nAt) by Vj". We notice that: if the random variable S” = Sj, then 
grt = S;41 with probability p 
S;-1 with probability 1 — p 


Therefore, the expected value of V at time step n at S; should satisfy 


vr = e TAt (p Vent (1 — p)vjrit) 


the expectation of the binomial distribution discounted by e~’4* 


option pricing is 


. Thus, the binomial method for 


(i) Choose N and define At = (T' — t)/N. N is chosen such that At < o?/r? 
(ii) Define Ax = ov At, and for each n = N,..., 0, define S; = Soe? with |j| < N. 


(iii) Define VN = A(S®), and define V;” forn = N —1, N — 2,...,0, by 


Vit =e (pvt + (1—p)Vi4}) , for |j| <n, j +m = even. 


(iv) Output: V, is the option price at time ¢ for asset price So. 
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Example. For put option, 


T = 5months = 0.4167 year 
At = 1months = 0.0833 year 
r = 01, c=04 
So = $50, EF = $50 
= e7V4t — 1.1294 
d = 0.8909 
D 0.5076 
et = 1.0084 


We begin to generate a binomial tree from Sp = 50 consisting of Si = Sou’d"—*, where n+j = 20, 
—n <j <n. Then we compute V;” inductively from n = N — 1 ton = 0 by 


Vit =e" (pyntt + (1 — pV). 


with ye being the payoff function A(S wv ). The value Ve is our answer. 
A C code for European option is given below. In computer code, we use the data structure 


e S[n|[m], V[n][m], m = 0,...,n, n = 0,...,N. Before, we use Si, where n + j = even and 
|j| < n. The relation between 7 and m ism = (j + n)/2. 


e The relation 


ie e Tat (pv! ser al — p)vjit) 


is replaced by 
Vin =e" (pVingi + (1 — p)Vin) - 


Binomial_European_option(V,S,S0,E,r,sigma,dt,N) 
{ 

double u,d,p; 

double discount, tmp; 

int m,n; 

discount = exp(-rxdt); 

u exp (Sigmax*saqrt (dt) ); 

d 1/u; 

Pp (1l/disecount — d)/(u-d); 


//Binomial tree for asset price: at level n. 
S[0] [0]=s0; 
for (n = 1; n <= N; ++n) 
{ 

for (m=n; m > 0; --m) 

S[n] [m] = uxS[n-1] [m-1]; 

S[n] [0] = d*S[n-1] [0]; 

} 


// Binomial tree for option. 
for (m=0; m <=N; ++m) 
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V[N] ae = payoff (S[N] [m],E); 
for (n= N-1; n >= 0; --n) 
{ 
for (m=0; m < n; ++m) 
Vin] [m] = discount«(p*V[n+1] [m+1] + (1-p)*V[n+1] [m]); 


5.3 Finite difference methods (for the modified B-S eq.) 


In this section, we shall solve the Black-Scholes equation by finite difference methods. Recall that 
the Black-Scholes equation is 


Ly g9PV _ av 


We set up the initial time is 0 and reserve the symbol ¢ for time parameter. Using dimensionless 
variables S = He*, V = Ev, 7 = T — t, we have 


1 ,07v il, Be 
t= 5° 5 + (r 5°) rv, 7T€ [0,7]. 
Let v = eu, then u satisfies 
_1 wo 1 9.08 
T ; € 0, T]. 5.6 
Ur = 507 + (r— 50%), rE [7 65.6) 


The initial condition for w is 


(20) = max{e* —1,0} for call option 
a max{1—e*,0} for put option 


The far field boundary condition for u is 
u(—oo,t)=0, u(x,t) = e%e"” as x 4 00 


for a call option, and 
u(—oo,t) = 1, u(x,t) =0 as % > oo 


for a put option. 


5.3.1 Discretization methods 


To solve numerically, we follow the following procedure: 


e Discretize space and time. We discretize [0,7] into N steps, Ar = T/N. We choose a 
small proper Az (to be determined later) and define x; = jAr,7 = —N,..., N. We shall 
approximate u(x;,nAr) by U?. 
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e Spatial discretization. We shall approximate the spatial derivatives by finite differences: 


— Uy is approximated by 
Ujt+1 — Uj-1 
2Ac 


Ux — 


— Ure is approximated by 
Ujt1 — 2u; + Uj-1 
(Ax)? 


Then the right-hand-side of is discretized into 


Urr * 


o Uyjay = 2U; + O51 Ltr o? Uae = Oy 


(QU); = (>) (Az)? P= oe 


e Temporal discretization. For the temporal discretization, we introduce the following three 
methods: 


— Forward Euler method: 


yet ye 
J I ny. 

we (QU). 

— Backward Euler method: 
prt _ n 

j _—_ n+1) . 

— Crank-Nicolson method: 
+1 
U. ; — ue 


ae = 5 llqu"*); + (QU). 


e Boundary conditions: we impose 
sae =0, Ut = eNAw rnAt 
for call option. 


For forward Euler method, we find U ; inductively in n forn = 1,...,N. We start from n = 0, 
which is the initial condition. We are given U}, j = —N,...,N. From nto n + 1, the boundary 
values idee and me are assigned. For —N < 7 < N, we find up by 


n+l. yrn ny). 


We stop this recursive process until all UN are found. 

For other methods, which involves Q(U"*') in the equation of U"+!. They can be expressed as 
a linear system for U"*!. There are standard method to solve such linear system. We shall discuss 
this issue later. 
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5.3.2 Binomial method is a forward Euler finite difference method 


We claim that the forward Euler method for wu in backward time 7 is indeed the binomial method for 
option price V if we choose Ax properly. Let us use V;” to approximate the solution V(a;,nAt), 
where At = Ar = T/N. We notice that the time variable for V is the forward time t and the 
time variable for u and v is the backward time 7. Thus, V(a;,nAt) = v(a;,(N — n)Ar) and 


v =e "Tu. We then have 
n _ (—r(N-n)Atr7N-n 
V; =e U; 4 
For the forward Euler method, we rewrite it in terms of V and also 


ee _ Ce ee cit Ar(QUu" **); 


In terms of eee we have 


rAty7n _ yrnt+l n+1 
styrene a AQ, 
1 At ge it 

—_ n+1 2 n+l1 n+1 n+1)\ , n+1 
VP Poet Vi 2p FE 


1 1 +1 
= aVjt + OVE + VE 


Ve) 


where 


At 
b = 2 n+l 
( (Ax)? Wes 
1 ANG. 6 = At 
c= r F 
2 (Ax)2” 2’ 2A 
Notice that 
a+b+c= 
We should require 
a,b,c > 0, 


for stability reason to be explained later. Thus V7” is the “average” of Vier ae ve with 
weights a, b, c, then discounted by eat 


The stability condition a, b,c > 0 reads 


2 
a | At AG 5 
<1 
| 2 | Az = (Az? = 
These give constraints on Ax and At: 
1 Ax)? 
Ax < Ae 
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Binomial case: b = 0 When we choose 
Az =oavAt, 


then b = 0. In this case, 
eae Ls pv ne (1 — p) Vpn 


where 
1 he ae At 
ae ae ee oe 
The stability condition is satisfied if and only if 
O0<p<l. (5.7) 


We see that this finite difference is identical to the binomial method in the previous section. 


5.3.3 Stability 


Definition 5.4. A finite difference method is called consistent to the corresponding P.D.E. if for any 
solution of the corresponding P.D.E., it satisfies 


ED.E (finite difference equation) + «(Azx, At) 


and « + Oas At, Ax > 0 


Definition 5.5. The truncation error of a finite difference method is defined to be the function 
e(Az, At) in the previous definition. 
For instance, the truncation error for central difference is 
2 2 
or o 
Ou= Zee +(r—- py ux + O((Az)?). 


And the truncation for various temporal discretizations are 


1. Forward Euler: 


u(jAz, (n+ 1)At) — u(jAz, nAt) 


a (Qu)(jAa, nAt) = O((Ax)*) + O(At). 


2. Backward Euler: 


u(jAg, (n+ 1)At) — u(jAz, nAt) 


(Qu) (jAz, (n + 1)At) = O((Ax)) + O(At). 


3. Crank-Nicolson method 


u(jAz, (n+ 1)At) —u(jAz,nAt) 1 
At 2 
= O((Az)’) + O((At)”). 


(Qu) GAa, (n+ 1I)At) + (Qu) GA, nAt)] 
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The true error Uj! — u(jAx, nAt) is usually estimated in terms of the truncation error. 


Definition 5.6. A finite difference equation is said to be (L? —)stable if the norm 


JU"? = So UPPAe 
j 


is bounded for all n > 0. 


Definition 5.7. A finite difference method for a P.D.E. is convergent if its solution U. a converges to 
the solution u(jAx, nAt) of the corresponding P.D.E.. 


Theorem 5.5 (Lax). For linear partial differential equations, a finite difference method is conver- 
gent if and only if it is consistent and stable. 


This theorem is standard and its proof can be found in most numerical analysis text book. We 
therefore omit it here. 

Since the consistency is easily to achieve, we shall focus on the stability issue. A standard 
method to analyze stability issue is the von Neumann stability analysis. It works for P.D.E. with 
constant coefficients. It also works “locally” and serves as a necessary condition for linear P.D.E. 
with variable coefficients and nonlinear P.D.E.. We describe his method below. 

We take Fourier transform of {U;}°° _.. by defining 


j=-CO 
O(E)= SD) Uje™ 
j=—-OCO 


It is a well-known fact that 


= A JOOP ae 
= | OI? 


bean 
j 


Thus, the boundedness of >), |U. ;|? can be estimated by using \|U7||?. The advantage of using U is 


that the finite difference operation becomes a multiplier in terms of U. Namely, 


Boe) = (Bee 


j 


3 — = oe) 
2 


| 
- (= ) s Uje8 
j 


(2isin €)U(€) 
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For the finite difference operator QU,m we have 


—_—.~ 


(QU) = J ~(Qu)je% 


= QU). 
For forward Euler method, 
Org) = (1+ AtQ(E)T 
= G()U" 
Gru 


We observe that 


[ \Fe@ra 


—pi 


/ G(OP"(e)2 a 


< _amax|G(6)P" / \T()/2 ae 


€€(—7,77) 


If |G(€)| < 1,)Vé € (—7,7), then stability condition holds. On the other hand, if |G(¢)| > 1 at 
some point €o, then by the continuity of G, we have that 


IG(@)| zlt+e 


for some small € > 0 and for all € with |€ — €| < 6 for some 6 > 0. Let consider an initial condition 


such that | | 
| 1). ee e6| 0 
0/¢) — 
Oey { 0 otherwise. 


Then the corresponding U” will have 


[ \Fe@ra = [ice rae? ae 


+> © 
as n —> oo. We conclude the above discussion by the following theorem. 


Theorem 5.6. For a finite difference equation with constant coefficients, suppose its fourier trans- 
form satisfies’ 


U"*1(£) = G(OT*(6) 


Then the finite difference equation is stable if and only if 


IG(g)| < 1 VE € (—7, 7]. 
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Example: Let apply the forward Euler method for the heat equation: uz = urz. Then 


perl _ ye 1 At 
Jj de! hoa n 9U” n n+l _ yn D? n 
At (Ane! ji — 2U; es Sw Oa: (Az)? e 
From von Neumann analysis: 
until = [1+ au 5 (2 cos € — 2)|0" 
(Ax)? 
BE Es 
= i 2>)U" 
ea (05 aa 
= Ge) 
Hence 4 i A ; 
t t 
IG(é)| <1 == (Any sin? : < —— (An? ae 5 (stability condition) 
If we rewrite the finite difference scheme by 
At At At 
(eh es Sey So 2= Ue Ue 
j (Ax)? jeit CL (ar gt (Ax)? j-l 


= OUF 4 OU a cu, 
Then the stability condition is equivalent to 
a,b,c> 0. 


Since we have a + b+ c = 1 from the definition, thus we see that the finite difference scheme is 
nothing but saying te is the average of U7, ,, U;’ and U;"_, with weights a, b,c. In particular, if 
we choose ar = 55 then b = 0. If we rename a = p, c = 1 — p, then Us = puri + (1- 

Pp) U;_,. This can be related to the random walk as the follows. 

Consider a particle move randomly on the grid points 7 Az. In one time step, the particle moves 
toward right with probability p and left with probability 1 — p. Let U;’ be the probability of the 
particle at 7 Az at time step n for a random walk. 

a = pu al = puss 

We can also apply the above stability analysis to backward Euler method and the Crank-Nicolson 
method. Let us only demonstrate the analysis for the heat equation. We left the analysis for the 
Black-Scholes equations as exercises. 


1. For backward Euler method, 


n+l n 
a 
At (AT 


(Up —2urt} + unt), 


5.3. FINITE DIFFERENCE METHODS (FOR THE MODIFIED B-S EQ.) 79 


Then K 
t —_ ae ——_ — 
(1 + (4sin? S) (pe = Ur) =) n th = G(€)U”, 
where ; 
G() = 
Lic 4 (sin? §) 


We find that |G(€)| < 1, for all €. Hence, the backward Euler method is always stable. 


2. For Crank-Nicolson method, 


De call ( 


Ae (Aa? (Ua — UF + Up) + Osea — 20; +U},)). 


jt+l1 


Its Fourier transform satisfies 


Ciara ee a2 &)\fjnti 2 S\ aa 
i = aa? | (4sin su (4sin 3U ; 
We have ‘K é x é 
| t 2 n+l t 2 n 
(14 thay Wiiiozs al hae (sin 5))U 
and hence 


Leta = 2 sin? s then G(€é) = 4£%. We find that for all a > 0, |G(€)| < 1, hence 


+a 


1 
Crank-Nicolson method is always stable for all At, Ax > 0. 


Exercise study the stability criterion for the modified Black-Scholes equation 


u- = an + (r aN 
T 9) 9 b) 


for the forward Euler method, back Euler method and Crank-Nicolson method. 


5.3.4 Convergence 


Let us study the convergence for finite difference schemes for the modified Black-Scholes equation. 
Let us take the forward Euler scheme as our example. The method below can also be applied to 
other scheme. 

The forward Euler scheme is given by: 


Ut — Ur 


og QU); 
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We have known that it has first-order truncation error, namely, suppose uy := u(jAx, nAt), where 
u is the solution of the modified Black-Scholes equation, then 


yntl _ yn 
We subtract the above two equations, and let e’’ denotes for w'—U; and ¢?' denotes for the truncation 
error. Then we obtain 


+1 
in er are 
aie (Qe"); + 6) 
Or equivalently, 
eet = aef 41 + be} + cep_1 + Ate} (5.8) 


Here, a,b,c > 0 anda+6+ c= 1. We can take Fourier transformation e” of e”. It satisfies 


ertl(€) = G(Q}E*(E) + Ater(€) 


where 
G(é) = ae® + b+ ce. 


Recall that the stability |G(€)| < 1 is equivalent to a,b,c > 0. Thus, by applying the above 
recursive formula, we obtain 


Jl < lent) + Aten 
< er? + At (|Ge>I| + lI) 
< er? + At (le + 41) 
< lle] + Ado Ne 
k=0 
< O(At) + O((Az)?). 


Here, we have used the estimate for the truncation error 
lle"|| = O(At) + O((Ax)?’, 
and that nAt = O(1). We conclude the error analysis as the following theorem. 


Theorem 5.7. The error er = u(jAz, nAt) — U. iM for the Euler method has the following conver- 
gence rate estimate: 


oy lepPAay < O(At) + O((Az)’), for all n. 
j 
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It is simpler to ontain the maximum norm estimate. Let E(n) := max; |e’’| be the maximum 
error. From (5.8), we have 


lee] < alessal + bel + clejal + Atle}| 
< aE(n) + bE(n) +cE(n) + Ate 
where 
e:= max|¢| = O(At) + O((Az)’). 
gn 
Hence, 


E(n+1) < E(n) + Ate. 


Since we take U. ; = ud, there is no error initially. Hence, we have 
n-1 

> Ate 

k=0 


nAte 


IA 


E(n) 


IA 


Since nAt is a fixed number, as we take the limit n — oo, we obtain the error is bounded by the 
truncation error. We summarize the above discussion as the following theorem. 


Theorem 5.8. The error e; = u(jAz, nAt) — U; for the Euler method has the following conver- 
gence rate estimate: 


max |e/| < O(At) + O((Az)?). 
a 
Exercise. Prove that the true error of the Crank-Nicolson scheme is O((At)”) + O((Az)?). 


5.3.5 Boundary condition 
For the modified Black-Scholes equation, we have 
u(—oo,t)=0, u(x,t) = e%e"” asx > 00 
for a call option, and 
u(—oo,t) = 1, u(x,t) = 0as x > co 


In computation, we can choose a finite domain (xy, xR) with rx, << —land ap >> 1. The 
boundary condition at the boundary points are an approximation to the above far field boundary 
condition. 

In practice, we don’t even use this boundary condition. Indeed, if we want to know u(xz, nAt), 
we can find the numerical domain of this quantity, which is the triangle 


{GAx,mAt) | |j—k| <n —m} 


We only need to compute u in this domain, which needs no boundary data. 
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5.4 Converting the B-S equation to finite domain 


The transformation « = log(S/E) converts the B-S equation to a heat equation. However, the 
domain of x is the whole real line. For numerical computation, it is desirable to have a finite 
computation domain. The transformation in this section converts S to € with € € (0,1). The 
price is that the resulting equation has variable coefficients. But this is not a problem for numerical 
computation. 

We define the transformation: 


S 
oS ae G2) 
— _ V(S,t) 
Ve= SLE (5.10) 
r = T-t. (5.11) 


Notice that € is dimensionless and important values of € are near 1/2. With this, the inverse trans- 
formation is 
Be. Ae eae 


MS Gasth ig = > 


We plug this transformation to the B-S equation. We allow o depend on S. Define a(€) = o( £€/1— 
€). Then the resulting equation is 


OV oil OV OV = 
Br = 3 (EA See + EL — Ge — r(L— 8)V, (5.12) 
for0 <€ < 1land7r > 0. The initial data reads 
sie gy 1 &y (HE 
V(E,0) = —z A (=) 2 (5.13) 


For a call option, the payoff is A(.S) = max(S — E70). The corresponding 


V(E,0) max(S — E,0)(1—£)/E 


max(2€ — 1,0). 


Similarly, V(€,0) = max(1 — 2€,0) for a put option. 

On the boundaries € = 0 and € = 1, the diffusion coefficients are degenerate. If the solution 
is smooth up to the boundaries, then on the boundary, the equation is degenerate to the following 
ordinary differential equations: 


OV (0,7) 
— = —rV (0, ie) 
OV (lise): 0 
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The corresponding solutions are 


V(0,r) = V(0,0)e"" (5.14) 
Vila) = Vial. (5.15) 


We can discretize equation (5.12) by finite difference method. Let A€ and Az are the spatial 
and temporal mesh sizes, respectively. Let €; = 7 A€, 7” = nAr. The boundaries points are €g and 
Ey. We use central difference for 0?V /O€? and OV /0E. The resulting finite difference equation 
reads 


dvj 1 o¢2 2 Vjt1 — 20; + Yj-1 
— = ~¢7¢*(1 
dr g 7585 (1 — &) Re? 
Uj — Uj 
bré;(1 — &)) A * — (1 - &)y; 


ORE 


We can discretize this equation in the time direction by forward Euler method. The stability con- 
straint is 
9| Ar 


2 
: <1. * 
: 1 (5.16) 


oAE = 


r&(1—§) — 5031 — &)?| So < FBO -§) 
Remark. Many options have non-smooth payoff functions. This causes low order accuracy for 
finite difference scheme. Fortunately, many simple payoff function has exact solution. For instance, 
the European call option. For general payoff function, we may subtract its non-smooth part for 
which an exact solution is available. The remainder is smooth, and a finite difference scheme can 
yield high-order accuracy. 


5.5 Fast algorithms for solving linear systems 


In the backward Euler method and the Crank-Nicolson method, we need to solve linear systems of 
the form 


AU =F. 
For the backward Euler scheme, 
A = diag (—a,1+a+c,—c) 
l+at+e —Cc 0 + 
—a l+a+ec —c 0 . 
0 —a l+a+ec -c 0 | 
0 —a 1l+a+e —Cc 0 
0 —a lta+ec —C 


84 CHAPTER 5. NUMERICAL METHODS 


and 
n n+1 
Ue +1 au; 
U= : = : 
n n+1 
eae cU;, 


For the Crank-Nicolson scheme We have 


AuU™t! a BU” an prti/2 


where A = diag (—§,1+ $+§,—§) B= diag ($,1— $—§,§), and 
a(UPt* + UP )/2 
0 


prtt/2 — 


OR 
+1 
Laue + Uf.) /2 
Now, we concentrate on solving the linear system 
Ax = f. 


The matrix A is tridiagonal and diagonally dominant. Let us rewrite A = diag (a,b,c). Here, the 
constants a,b,c are different from the average weights we had before. We may assume b > 0. 
We say that A is diagonally dominant if b > |a| + |c|. More generally, A may takes the form 
A = diag (a;,b;,c;). and |b;| > |a| + |c;|. Without loss of generality, we may normalize the 
j — th so that b; = 1. 

There are two classes of methods to solve the above linear systems. One is called direct methods, 
the other is called iterative methods. For one-dimensional case as we have here, direct method is 
usually better. However, for high-dimensional cases, iterative methods are better. 


5.5.1 Direct methods 
Gaussian elimination 


Let us illustrate this method by the simple example: A = diag (a,1,c). We multiple the first 
equation by —a and add it into the second equation to eliminate the term x;,+) in the second 
equation. Then the resulting equation becomes 


1 Cc 0 0 .:-:- 0 Vir tl bj, 41 
0 l-ac c 0 ::: O Lip +2 —abj,41 + bj, +2 
0 a 1 c¢ 0 Vip 43) = 5.43 

0 : 


1} |&jp-1 Djp-1 
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We continue to eliminate the term a in the third equation, and so on. Finally, we arrive 


1 ¢ 0 0 O] fayjp4i bad 
0 1l-ac @ OL see 0) |e bi oe 
0 0 1—c/(l—ac) c O} | 243] = | 55,43 


0 0 0 : 4G) : : 
0 0 0 Os Sette (hi) eee g bi 1 


Then x; can be solved easily. The diagonal dominance condition guarantee that the reduced matrix 
is also diagonally dominant. Thus, this scheme is numerical stable. 


LU decomposition 


We decompose A = LU, where 


[Or ~Oe. eey 10 okie Wee We SO 
a 1”. : | 0 Ujp+2 7 : 
L= | a a: rs 0 = | GeO, | 
, a | 5. Pit | 

ie Orta 0: sO) 0 ayes 


It is easy to find a recursion formula to find the coefficients @, u and v’s. Once these are found, we 
can find x by solving 
Ly =6b,Ux = y. 


These two equations are easy to solve. One can show that both L and U are diagonally dominant if 
Ais. 
If we watch carefully, LU-decomposition is equivalent to the Gaussian elimination. 


Cyclic reduction method 


Let us take the case A = diag (a, 1, c) to illustrate this method. Consider three consecutive equa- 
tions 


ON ig 9 Fon 4 baiog -S- banat 


GXq;-1 + Lo; + c®oj41 = 09; 


AX2; + Laj41 + CLaj+2 boj44 


We can eliminate the odd-index terms x2;-; and x2;41. Namely, —a x (2j — 1)—eq +(2j)-eq 
—c x (27 + 1)-eq: After normalization, we obtain 


! ! ! 
QXQj-2 + Laj + C LQj42 = b; 
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Here, 


2 2 
n a , c 
a= 


7 Cc se ’ 
1 — 2ac 1 — 2ac 
bi, = (ba; — abaj—1 — cbaj41)/(1 — 2ac). 


If we rename 2’, = x9; Then we have A’x’ = b’, where A’ = diag (a’,1,c’). Notice that the 
system is reduced to half and with the same form. One can show that the iterative mapping 


A 4 

rH], 

c c 

converges to (0,0)! quadratically fast, provided |a| + |c| < 1 initially. Thus, for few iteration, the 


matrix A is almost an identity matrix. We can invert it trivially. Once x2; are found, the odd-index 
x + 23 + 1 can be found from the equation: 


QE; + Laj41 + c&a;49 = bo541. 


A careful reader should find that the cyclic reduction is also a version of the Gaussian elimination 
method. 


5.5.2 Iterative methods 


Most iterative methods can be viewed as a proper decomposition of A, then solve an important and 
treat the rest as a perturbation term. 


Jacob method 


In Jacobi method, wer decompose 
A=D+B 


where D is the diagonal part and B is the off diagonal part. Since A is diagonally dominant, we 
may approximate x by the sequence x”, where x” is defined by the following iteration scheme: 


Da"t! + Br” = b. 
Let the error e” := x”"*! — x”. Then 
De” = —Be”"! 


Or 
e"=—D'Be™!. 


Let us define the maximum norm 
len | = es le; 
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Then 
_ @n-1 __ © n-1 
lef] = |- 5-1 a pcitl 
lal may y Il on 
< Lyle “le jal ‘| 
lal + lel on 
a lel 
Hence, 
lle" < alle" | 
where p = allel < 1 from the fact that A is diagonally dominant. This yields the convergence of 


[6] 
the sequence x,,. The limit x satisfies the equation Ax = b. 
Gauss-Seidel method 


In Gauss-Seidel method, A is decomposed into A = (D + L)+U, where D is the diagonal part, L, 
the lower triangular part, and U, the upper triangular part of A. The approximate solution sequence 
is given by 

(Baha +o" = 6. 


As before, the error e” := x"t! — x” satisfies 
n —1 n-1 
e” = -(D+L) ‘Ue 


To analyze the decay of e”, we use Fourier method. Let 
en(é) := S- ere 8, 
j 


Then we have 


ee) = Geyer, 
ces 
b+ ae 


It is easy to see that the amplification matrix G' satisfies 


max |G(é)| := p <1, provided |b| > |a| + |c]. (5.17) 


This shows that the Gauss-Seidel method also converges for diagonally dominant matrix. 
Exercise. Show the above statement (5.17). 
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Successive over-relaxation method (SOR) 


In the methods of Jacobi and Gauss-Seidel, the approximate sequence x” is usually convergent 
monotonely. We therefore have a chance to speed them up by an extrapolation procedure described 
below. 


yore —D7"(L+U)x" +b 


gh = gt a og _ a”), 


Here, w is a parameter. In order to speed up, we require w > 1. We also need to require w < 2 for 
stability. The optimal w is chosen to minimize the amplification matrix G,,(€). 
Exercise. Find the amplification matrix G,, and the optimal w for the matrix A = diag (a, b,c). 
Also, determine the rate 
p := min max |G,,(§)]. 
wg 


Chapter 6 


American Option 


6.1 Introduction 


An American option has the right to exercise any time during the life of the option. We denote the 
American call by C and American put by P. The first important thing we should note is that the 
value of an American option is greater than or equal to the payoff function: V(S,t) > A(S,t). 
Otherwise, there is an arbitrage opportunity because we can buy the American option then exercise 
it immediately to gain a net profit V — A. 

We recall that the value of an American call option is equal to that of a European call option, 
see Theorem 3.2} However, for other cases like the the American put option or the American call 
option on dividend-paying asset, the American options do cost more. We explain why it is so below. 


American put First, we notice that there must be some region in 0 < S < EF in the S-P plane 
where the European put p(S,t) < E— S. In fact, we know that p(0,t) = Ee-"?-9 < E-0, 
which is below the line p = E — S ina neighborhood of S' = 0, see Figure[6. 1] On the other hand, 
the corresponding American option P(S,¢) must satisfy 


P(S,t) > A(S)=E-S. 


Otherwise, we can buy an asset S and a put P(S,t), then exercise it immediately. We have the 
right to sell S on E. This gives a risk-free profit: # — P —S > 0. Thus, in the region where 
p(S,t) < E—S < P(S,t), we must have p(S,t) < P(S,t). 


American call on a dividend-paying asset For an European call on dividend-paying asset, its 
value must lie below the payoff function A(.S) = S — F in some region S > E, see Figure 6.2. In 
fact, c(S,t) ~ Se Pot -*) for large S, which is below the line c = S — E. Therefore, there must 
be a region in S > E in (S,t)-plane where c(S,t) < A(S) := max{S — E,0}. In this region, if 
we could exercise the corresponding American call option, then based on V(S,t) > A(S,t) (ie. 
C(S,t) > S — E), the corresponding American call option C' must also satisfy 


C(S,t) > A(S) = max{S — E,0} > c(S,t). 
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Ee"(T-) 


Figure 6.1: The dash line is the price of an European put option. There is a region where p(S,t) < 
max{E — S$,0}. This can be seen from the boundary condition P(0,t) = Ee~"?—, which is 
below the payoff function p = E — S. 


We summarize the discussion as below. 
Proposition 6. American option V always satisfies V(S,t) > A(S). Furthermore, 
e for American put option, P(S,t) > p(S,t) in some region where p(S,t) < Ap(S); 
e for American call option on dividend-paying asset, C(S,t) > c(S,t) in some region where 
NGS), 
6.2 American options as a free boundary value problem 


6.2.1 American put option 


Existence of optimal exercise price S;(t). First, there must be some value of S for which it is 
optimal from the holder’s point of view to exercise the American option. Otherwise, we should 
hold the option for all possible S. Then this option is identical to a European option. But we have 
seen that this is not the case. Hence, for each time ¢, there is an optimal exercise price S(t) such 
that, when S < S'r(t) one should exercise the put option, and when S > S'r(t), we should hold the 
option. In other words, we should have: 


e P(S,t) = max(H — S,0) for S < S(t), 


e P(S,t) satisfies the Black-Scholes equation for 5 > S(t). 
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Figure 6.2: Dash line is the European call option on a dividend-paying asset. Notice that c(S,t) ~ 
Se~PolT-9 for § >> 1, hence C(S,t) < A(S) = max(S — E,0) for some S'f(t). 


The first statement can be viewed as the definition 5';(t). That is, 
Sy(t) := max{S|S < EF and P(S,t) = max(E£ — S,0)}. (6.1) 


The optimal exercise price can also be viewed as the minimal exercise price that one should hold 
the put option: 


S(-) = min{b(-)|P(S, t) is a B-S solution in S > b(t) satisfying P(b(t),t) = E — b(t)}. 


This means that, given a exercise price curve b(-), the holder holds the option if S > b(t) and 
exercise it if S < b(t). When it is exercised, the payoff is E — b(t). The value of P(b(t),t) must 
be equal to FE — b(t) by the no-arbitrage assumption, otherwise we can buy or short P then exercise 
it right away to gain a net profit. Thus, given an exercise price curve b(-), we can find a solution 
P(S,t) with boundary condition P(b(t),t) = E — b(t). The optimal exercise price S,(-) is the 
smallest such exercise price curve. 
Notice that 
S$ F(t) <E. 


If S (or S¢(t)) goes up beyond £, then the right of a put to sell a higher price S with a lower price 
E is meaningless. Thus, such put P(.S,t) has no value at all. This means that the optimal exercise 
price cannot exceed FE. Furthermore, at time T’, we should have 


S;(T) = E. 


This follows from the Definition (6.1). 
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However, we don’t know the optimal exercise price S(t), t < T in advance. It is a new 
unknown, called free boundary. On this free boundary, there will be two boundary conditions for 
P. One boundary condition is used to determine the B-S solution Pgs uniquely in the region 
S>S f (t),0 <t < T. The other boundary condition is used to determine the free boundary itself. 
The following proposition gives the proper boundary conditions on S'¢(t). 


Proposition 7. Based on the no-arbitrage hypothesis, the price of an American put should satisfy 
the following boundary conditions on the free boundary S's (t): 


e P(S;(t),t) = BE — Sy(t), 

e SF (Sp(t),t) = —1. 
Remark. Notice that in the region S < S,(t), P(S,t) = E — S, where we have 0P/05(S,t) = 
—1. Therefore, the above boundary conditions are equivalent to 

e P(S,t) is continuous across the free boundary S = Sp(t). 

e oF (S, t) is continuous across the free boundary S = S(t). 


Proof. 1. Given the optimal price S'f(t), it means that it is the price to the holder to excise it, 
the return is £ — S(t). The corresponding option price must equal EL — S¥(t). Otherwise, 
there is an arbitrage opportunity. Thus, we get P(S(t),t) = EB — S;(t) forO <t <T. 


2. The free boundary S = S;(t) is an unknown. We want to derive a constraint to determine it. 
Suppose we are given a boundary S = b(t), 0 < t < T. We solve the B-S equation in the 
region S > b(t), 0 <t < T with boundary condition 


P(b(t), t) = E — v(t) 


and final condition P(S,7T) = max(£ — S,0). Such a solution depends on the path b(-). We 
denote it by P(S,t; b(-)). The optimal exercise price S(t) is the minimal such curve D(-). 
Below S(t), then one should exercise the put option. Therefore, a small variation of P with 
respective to the change of b(-) is zero at the optimal exercise price S'r(-). This means 


5P(S/(t),t35¢(-)) 
5b(:) 


= 0. 


Here, 6P/5b means the variation of P w.r.t. the path b(-). 


'The variation of P w.r.t. a path b(-) is defined as below. Consider a small change of the path, say b(-) + er(-). The 
function r is called the direction of the variation. Then, just like the definition of gradient in n dimensions, we define the 
variation of P w.r.t. b(-) at a path b(-) by 


< |e=o P(S,t;b + er) = Se ds. 
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3. On the boundary S = b(t), P satisfies the boundary condition 
P(b(t), t; b(-)) = E — b(t). 


We can take variation of this formula w.rt. b(-), then evaluate at the optimal exercise price 
S(-). This gives 


OP 5(. 5P(S;(t),t; Sp(- 
55 SHO SO) SY b-)=5;() 4 \ oe HD 5 


Since the second term on the left-hand side is zero. We get 


PF (Sp(0). 6870) = 1 


Thus, we treat the American put option as the following free boundary value problem. There exists 
an optimal exercise price S'f(t) such that 


1. for S < S¥(t), early exercise is optimal, and P(S,t) = E—S; 
2. for S > S(t), one should hold the put option and P satisfies the Black-Scholes equation: 
OP... 1 OP OP 
Lg 6? 


aE + a? © age ag TF =o 


3. across the free boundary (S(t), t), both P and ge are continuous. 


We will solve this free boundary value problem in the next section. 


6.2.2 American call option on a dividend-paying asset 


As we have seen in the introduction of this chapter that an American call option C(S,t) on a 
dividend-paying asset has asymptotic value C(S,t) ~ Se~?0(7—*) for large S'. This value is below 
the payoff function A = max(S — E,0). Therefore, there must an optimal S(t) such that we 
should exercise this call option when S > S(t) and hold it when S < S¥(t). On the free boundary 
S' = S;(t), based on the no-arbitragy hypothesis, we should require that 


e both C and 0C/OS are continuous across the free boundary (S(t),t) forO << t < T. 
We summarize this by the following equations 
O° 50°C OC 
Ss? | Do)S 
z° age t — PulS ae 
C(S,t) = A(S) = max{S — E,0}, S > S(t). 
On the free boundary S = S'-(t), the boundary condition is required 
C(S(t), t) = Sp(t) — B, 


C4 ro = 0,0 48 S50) 


OC 
sas 1. Veer 
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6.3 American option as a linear complementary problem 


The American option can also be formulated as a linear complementary problem, where the free 
boundary is treated implicitly and numerical schemes can be easily designed. To illustrate this linear 
complementary problem, first we notice that an American option V should satisfy the following 
conditions: 


(i) V(S,t) > ACS), 


Gi) Yi + 4028? SY <r(V— = 


(iii) either V(S,t) = A(S) or V, + 4079? ov =r(V—- SX) should hold, 
(iv) both V and oy are continuous. 


Here, A(.S') = max(E — S,0) is the corresponding payoff function. 

We have seen the reasons for (i), and (iv). We explain the reasons of (ii) below. Let us consider 
the portfolio I] = V — AS. As we have seen that the Delta hedge eliminate the randomness of II 
and yields 


dadl=Vz+ a” S Dee" 
When it is optimal to hold the option, then 
OV 
dil=r(V —- = S 
r(V — 555). 
Otherwise, we should have i 
V 
dil < r(V — —S 
r(V — 555), 
based on no arbitrage opportunities. Thus, the Black-Scholes is replaced by the Black-Scholes 


inequality. 

To show (iii), we know that if we exercise the option, then V = A, otherwise, we hold the option 
and its value should satisfy the Black-Scholes equation. 

Properties (i)-(iii) can be formulated as the following linear complementary equations: 


V-A>0 
Vi + 4079? S¥ (Vv - Sa5) <0, (6.2) 
(V — A) (vi+4 1 1g2g22V _ wy s3t)) =i 


332 


We look for C! solution (i.e. both V and ov are continuous) in the domain S € (0,00) andt > 0 


with final condition 
PCS. T= ACS): 


Such a problem is called a linear complementary problem. The advantage of this formulation is that 
the free boundary is treated implicitly, which can be solved by some penalty method or projection 
method to be explained below. 
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Remark For American put, in the region S < S(t), P(S,t) = E —S. We plug it into the 
Black-Schole equation 


Thus, the B-S inequality holds. 


6.4 *Penalty method for linear complementary problem 


We can reformulate this problem in terms of x variable. As before, we use the following change of 
variables: V = Ev, S = Ee”, 7 = T —t. The free boundary now in «-variable is x(t). The free 
boundary value problem is formulated as 


(i) for —oco < & < xf(t), v = 1 — e*, and 


2 2 
o o 
Ur 5 Vex — (rT 5 juz t+ rv > 0. 
(ii) for x(t) < x < co, u > 1 — e”, and 
2 2 
Ur = Vex — (r = juz trv = 0, 


(iii) both v and gu are continuous. 
The linear complementary problem is formulated as: 
(i) v-(1-e*) 50, 


o2 


(ii) vu; gq Vara (r 5} Ux Prue ws; 


(iii) (vp — Sven — (r — vy + rv)(v — (1—e*)) =0, 


(iv) v and vz are continuous. 


with initial condition v(7,0) = A(z) =1-e”. 
We may replace v by ue~"” to eliminate the term rv. Then we have 


2 2 
oO oO 
(u 5 Ure — (7 5 Jue) (u-g) = 0 
2 2 
o o 
Ur 5 Use — (Tr 5 juz > 0 
U— g > 0, 


pa 
8 
© 
| 
S 
8 
= 
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with u, uz; being continuous. Here g(x,7) = max{e’’(1 — e”), 0}. The far field boundary condi- 
tions are 
u(x, 7) — 0, as Z > 0, Ulrer) rE er as ZT —> —OOo. 


A mathematical theory called parabolic variational inequality gives construction, existence, unique- 
ness of the solution. (see reference: A. Friedman, Variational Inequality). Let us demonstrate this 
theory briefly. The method we shall use is called the penalty method. Let us consider the following 
penalty function: 

on(v) = —eN”, 


It has the properties: (i) éy > 0, (ii) év(v) — O whenever v > 0.. We consider the following 
penalized P.D.E.: 


jue +on(u-g) = 0 
u(x,.0). = o(#;0), 
with u + 0asx—--+o00o, u—e’ asx —- —oo. 
From a standard theory of nonlinear P.D.E. (by monotone method, for instance), one can show 


that the solution uy exists for all NV > 0. We then need an estimate for uy and SUN The bounded- 
ness of these two gives that uy has a convergent subsequence, say uy, such that 


UN; — U, 
with wy, U, Zu N;> Ux being continuous. Moreover, 
|dn, (un, — 9)| < constant 
As N; — oo, we conclude u — g > 0. Further, on the set {u — g > 0}, On, (un, — g) + 0. Hence 
we have 


1 o? 
Ur — 5 Ulex (r 5 uc = Oon {u—g > O}. 


6.5 Numerical Methods 


6.5.1 Binomial method for American puts 


We recall that the solution of the Black-Scholes equation can be discretized by the following bino- 
mial method (or the forward Euler method): 


e Choose a large N, define At = T/N, Ax = oV At, Ss} = Spei*, lj] <n, 7 +n = even. 


e Define 
iad i = DV 4 wo’ 
where 
= ~ ya eA da 1/u,qg=1-—p. 
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Similarly, the linear complementary problem can be discretized as 
At 1 n+1 

(evs — pvp! — avin) (YP— AS) = 0. 

cay = BVA _ av) 

ya 


IV IV 


Here, A? := E — S®. This recursive procedure finding V" from Vv"! forn = N—1,N —2,...,0 
with initial condition 
N _ AN 
V5 — me . 
This discretized linear complementary problem can be solved by the following projected forward 
Euler method: 


Vit = max{e"™™ (pV) + Vj"), AFF. ce 


Here, A? := A(S®). This is called projected binomial method. 


Theorem 6.9. The projected binomial method for American option is stable and converges. Namely, 
VP > V(S,t) with Sj = SoeI4* + S and nAt — t as n —> oo. Furthermore, V(S,t) satisfies 
the linear complementary equations (6.2). 


Sketch the idea of the proof. 


1. One can show that this method is stable and converges in C} (IR*), the space of C!-continuous 
and bounded functions, by using the monotone operator theory. (Ref. Majda & Crandell, 
Math. Comp. ) An operator ® : Cy(IR) — C,(R) is called monotone operator if f > g 
implies ®(f) > @(g). For continuous case, the solution operator ®(t) of the Black-Scholes 
equation is a monotone operator. Their discrete version, the one-step solution operator which 
maps V"*! +5 V™” is also a monotone operator. The projection operator is a monotone 
operator, which can be proven easily by the fact that max(a, c) < max(b,c) ifa < b. 


2. Furthermore, from the projection V;" = maxf{e"*(pvrat + qVj4'), Avy, we obtain 


j+1 
pea 
—rAt +1 +1 

Vee? (ov + Vir) ) 
Thus the limiting function V(S,t) also satisfies V(S,t) > A(S,t) and the Black-Scholes 
inequality. 

3. From construction, the projection: V;" = max{e" (pvt + qv"), A’}} satisfies either 

V= A} = 0, or Vi" — ett (Vat + qv") = 0. Thus, it satisfies the complementary 
condition: 


rAtyyn n+1 n+1 n ny) _ 
(e*vp — pvp? — ave) (YP — a9) =0. 


From convergence result, their limiting function also satisfies the continuous complementary 
condition. 
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Below is a C code for binomial method for American put. 


Binomial_American_option(V,S,S0,E,r,sigma, dt,N) 
{ 

double u,d,p; 

double discount, tmp; 

int m,n; 

discount = exp(-rxdt); 

u = exp(sigmaxsqrt (dt) ); 

d = 1/u; 

p = (1/discount - d)/(u-d); 


//Binomial tree for asset price: at level n. 
S[0] [0]=S0; 
for (n = 1; 
{ 


n <= N; ++n) 
for (m=n; m > 0; --m) 
S[n] [m] = u*S[n-1] [m-1]; 
S[n] [0] = d*S[n-1] [0]; 
} 


// Binomial tree for option. 
for (m=0; m <=N; ++m) 
VIN] [m] = payoff_put(S[N] [m],E); 
for (n = N-1; n >= 0; -~-n) 
{ 
for (m=0; m < n; ++m) 
{ 
tmp = px*xV[nt+1] [m+1] + 
tmp *= discount; 
V[{n] [m] = max(tmp, payoff_put (S[n] [m],E)); 


(1-p) *«V[n+1] [m]; 


} 


double payoff_put (S,E) 
{ 
return max(E-S,0); 


} 


6.5.2 Binomial method for American call on dividend-paying asset 


The linear complementary problem for this American call option is 


(i) V(S,t) > A(S) = max{S — E,0} 


Gi) y+ 5 S2¥ + (r 


Do) S83 —rV <0, 


(iii) (Vi pe ree er DiS rv) (V—A)= 


(iv) V and Vg are continuous. 


The binomial approximation for the B-S equation is 


0. 


rAty7n __ n+1 n+1 
eV = iVj41 + avin, 
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where 


ao” At ao? At ao” At G7 
DS oqAny ee poke age Oe 


and p + q = 1. We choose At and Az so that p > 0 and q > 0. For American option, V has to be 
greater than A(S,t), the payoff function at time t. Hence, we should require 


Vj = max{e "(pV + qVj"4"), AF}. 


The above is the projective forward Euler method. 
For the corresponding binomial model, first we determine the up/down ratios wu and d for the 
riskless asset price by 


pu + (1 _ p)d = er oat pu a (1 = p)d? = e(2r—Do)+o?)At 
or equivalently, 


u = A+VA2-1, d=1/u, 


1 

A = sore ty el Doro? )Aty 
e(r—Do)At d : 

P u— d ’ qd 7 p 


Then the binomial model is given by 


gn use, with probability p 
i \ dst, — with probability q 
Ss, =s 


and 


Vp = max{e™™ (pV + aVj"4'), Sp — EY. 


6.5.3 *Implicit method 


For implicit method like backward Euler or Crank-Nicolson method, we need to add the constraints 
u®t! > g”*1 for American option. It is important to know that if an iterative method is used, then 
we should require this condition hold in each iteration steps. For instance, in the SOR iteration 
method, 


yr(ktl) max{V™() + ww (yt) = az) ars th 
es DD) Ue afer ak 0 eile 
Ve a Vee 


This guarantees that V”+! > Ar+1, 
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6.6 Converting American option to a fixed domain problem 


6.6.1 American call option with dividend paying asset 
We consider the American call option on a dividend paying asset: 
GOV OV 
S$? Do)S 
2° age tr Poids 
V(S,t) > A(S) = max{S — E, 0}, 
Cao ASA Ot ae 


V; 4 rV =0, 


where S(t) is the free boundary. On this free boundary, the boundary condition is required 
V(S4(t),t) = S4(t) — B, 


OV 
—— = <t<T. 
ag (Sit)st) =1, OStST 


We also need a condition for S f at final time 
S»(T) = max(L,rE/Do) 


Before converting the problem, we first remove the singularity of the final data (i.e. non-smoothness 
of the payoff function) as the follows. We may substract V by an European call option c with the 
same payoff data. Notice that c(S,t) has exact solution. The new variable V — c satisfies the same 
equation, yet it has smooth final data. 

To convert the free boundary problem to a fixed domain problem, we introduce the following 
change-of-variables: 


g = S/Ss(t) 

7 = T-t 

u(,7) = (V(S,t) —e(S,t))/E 
s(t) = S,(8)/E 


The new equations for these new variables are 


Be = FORE + (7 Do) + GL) eM ru, OSES, 


s¢ dr 
0<7T<T, 
u(l,T) = 9(s#(7),7); 0<7r<T 
Se (1,7) = A(s;(r),7), O<7<T 
s (0) = max(1,r/Do), 
where 
gls¢(t)y7). =. spt) — 1 — cE s7(7),T —7); 


h(sg(r),7) = sg(r)|1— SS (Bag(r),7 7) 
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At the boundary € = 0, the Black-Sholes equation is degenerate to 


Ou 
— =-TUu. 
OT 
With the trivial initial condition yields 
u(0,7) =0,0<7<T. 
In practice, we can solve the modified Black-Sholes equation with the boundary conditions 


{ u(0,7) =0 0<7<T (65) 


Be (1,7) =A(ss(7),7), OS TST 


We can differentiate the other boundary condition in 7 and yield an ODE for the free boundary: 


with s (0) = max(1,r/Do). 


6.6.2 American put option 


For American put option P, the B-S equation is on the infinite domain S > S f(t), 0O<t<T. 
Through the change-of-variable 


ae EP(S. 
t 
u(n,t) = 2°60 


the infinity domain problem is converted to a finite domain problem: 


=~ 2 


0<t<T 
u(n,T) =max(y7—£,0), O<n<n,(T), 
ulns(t), )=ngs(t) — £, 02227, 
an (np (t), t) = 1, O<t<T. 


We can further convert it to a fixed domain problem as that in the last section. 
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Chapter 7 


Exotic Options 


Option with more complicated payoff then the standard European or American calls and puts are 
called exotic options. They are usually traded over the counter. Their prices are usually not quoted 
on an exchange. We list some common exotic options below. 


1. Binary options 

2. compound options 
3. chooser options 
4. barrier options 

5. Asian options 

6. Lookback options 


In the last two, the payoff depends on the history of the asset prices, for instance, the averages, 
the maximum, etc., we shall call these kinds of options, the path-dependent options, and will be 
discussed in the next Chapter. 


7.1 Binaries 


The payoff function A(S) is an arbitrary function. One particular binary option is the cash-or- 
nothing call, whose payoff is 
A(S) = BH(S — E). 


This option can be interpreted as a simple bet on an asset price: if S > F at expiry the payoff is B, 
otherwise zero. We have seen its value is 


V ae | P(S',T, S,t)A(S')dS! = e "7-9 BN (da), 
0 
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where ; , 
1 log(S)— (r= 5 )(T=t) 
PS. 84) — 202(T—t) 
2no?(T —t) 


is the transition probability density for asset price in risk-neutral world. 


7.2 Compounds 


A compound option may be described as an option on an option. We consider the case where 
the underlying option is a vanilla put or call and the compound option is vanilla put or call on 
the underlying option. The extension to more complicated option on more complicated option is 
relatively straightforward. There are four different classes of basic compound options: 


1. call-on-call, 

2. call-on-put, 

3. put-on-call, 

4. put-on-put. 
Let us investigate the case call-on-call. Other cases can be treated similarly. The underlying option 
: Expiry: To, Strike price: Eo. 
The compound option on this option 

Expiry: T; < J, Strike price: F}. 


The underlying option has value C'(S, t, T2, Ez). At time T}, its value C(S, 7), T>, Ea). The payoff 
for the compound call option is max{C(S, 7,72, £2) — E1,0}. Because the compound options 
value is governed only by the randomness of 5, according to the Black-Scholes analysis, it also must 
satisfy the same Black-Scholes equation. We then solve the Black-Scholes equation with payoff 


max{C(S, 7), T2, E2) — E41, 0}. 


7.3. Chooser options 


A regular chooser option gives its owner the right to purchase, for an amount F at time 7}, either 
a call or a put with exercise price EF» at time T. Thus, it is a “call on a call or put”. Certainly, we 
have T; < T>. The payoff at T; for this call-on-a-call-or-put” is 


A= max{C(S, T\) — Fi, P(S, T) = Fy, O}. 


The compound option also satisfies the Black-Scholes equation for the same reason as above. From 
this and payoff function at T,, we can value V at t. The contract can be made more general by 
having the underlying call and put with different exercise prices and expiry dates, or by allowing 
the right to sell the vanilla put or call. By using the Black-Scholes formula for vanilla option, there 
is no difficulty to value these complex chooser options. 
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7.4 Barrier option 


Barrier options differ from vanilla options in that part of the option contract is triggered if the asset 
price hits some barrier, S = X, say at some time prior to T. As well as being either calls or puts, 
barrier options are categorized as follows. 


1. up-and-in: the option expires worthless unless S reaches X from below before expiry. 
2. down-and-in: the option expires worthless unless S reaches X from above before expiry. 
3. up-and-out: the option expires worthless if S reaches X from below before expiry. 


4. down-and-out: the option expires worthless if S reaches X from above before expiry. 


7.4.1 down-and-out call(knockout) 


A European option whose value becomes zero if S ever goes as low as S = X. Sometimes in the 
knockout options, one can have boundary of time, or one can have rebate if the barrier is crossed. 
In the latter case, the option holder receives a specific amount Z for compensation. 

Let us consider the case of a European style down-and-out option without relate. We assume 
X < E. The boundary conditions are 


V(X,t) =0, (boundary condition), V(S,t) ~ Sas S — oo. 


The final condition, V(.S, 7) = max{S — E,0}. For S > X, the option becomes a vanilla call, it 
satisfies the Black-Scholes equation. 
Let us find its explicit solution. Let S = Ee*,t = T — GIDE V = Ev. The Black-Scholes 
equation is transformed into 
Ur = Vee + (Kk — 1)vz — kv, 


where k = ss We make another change of variable : 


y = ett Pry, 


We choose a, ( to eliminate the lower order terms in the derivatives of x: 


BerttPTy 4 eH PTy = ate tPTy + Qe tOry, + OPH PT a 
+(k — 1)(ae™ try + et P Ty) — ket Pry, 


This implies that a = —3(k — 1) and 8 = —}(k + 1)? and equation becomes u; = Une. 
Let 2p = log(4), or X = Ee*°. The boundary condition becomes u(x%o,7) = 0 and u(x, 7) ~ 


e(l-a)a—Br ag x 5 oo. The initial condition becomes 
ul(e;0) ule) = max{e2 (hth) — e2(k-1e o1. 


This follows from the payoff function being 


A =max{S — E,0} = Bmax{2 — 1,0} = Emax{e” — 1,0}, 
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and V = ewe tPru(e, Trea); with u(x,0) = uo(x) = e~* max{e*—1,0} = max{el 


e °” 0}. Notice that because X < E, we have xo < 0, and 


—atl)x_ 


(x) 0 forzgp <x <0, 
MOT V max{e(-e+De — e-@ 9} otherwise. 


We use method-of-reflection to solve above heat equation with zero boundary condition. We reflect 
the initial condition about xo as 


_ f uo(z), for %p < © < 00 
u(x, 0) = { —ug(2%9 —2), for —co< a4 < Zo. 


The equation and the initial condition are unchanged under the change-of-variable: x — 22%9 — x, 
u — —u. From the uniqueness of the solution, the solution has the property: 


u(2x9 — x,t) = —u(z, t). 


From this, we can obtain that u(xo,t) = —u(ao,t) = 0. 
Since C = Ee®*+?7 yu, is the vanilla call, where wu satisfies the heat equation with the initial 


condition: 


k+1)a k—-1)ax 


for x > 0 
0, forz <0 


1 1 
2 _ e2( 


U1 (x, 0) = { 
Using this and the method of reflection, we may express V in terms of C as the follows. First, we 
may write V = & enter (u; + ug), where the initial condition for uz is the reflected condition from 
U1: 
_ f uo(2%9-— 2,0) forx <0 
ual, 0) = { 0 for xz > 0 
The solution u ; corresponds to C(S,t). The solution ug is corresponds to 
e20(%—20) C(x? /S, t). We conclude 


V =C(S,t) - (Fy Y0(x?/S, 0) 


7.4.2 down-and-in(knock-in) option 


An “in” option becomes worthless unless the asset price reaches the barrier before expiry. If S 
crosses the line S = X at some time prior to expiry, then the option becomes a vanilla option. It is 
common for in-type barrier option to give a rebate, usually a fixed amount, if the barrier is not hit. 
This compensates the holder for the loss of the option. 

The boundary condition for an “in” option is the follows. The option is worthless as S' — oo, 
ie., V(S,t) > 0as S > 00. AtT, if S > X, then V(S,T) = 0. Fort < T, V(X, t) = C(X,t). 
Since the option immediately turns into a vanilla call and must have the same value of this vanilla 
call. For S < X,V(S,t) = C(X, t). We only need to solve V for S > X. V still satisfies the same 
Black-Scholes equation for all 5’, t, because its randomness is fully correlated to the randomness of 


S. 
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We may write V = c — V, where c is the value of a vanilla call. Then the boundary condition 
for V is V(S,t) =c-V~S-—0=SasS — oo. And 


V(X,t) =c(X,t)-V(X,t) =0, V(S,T) =e(S,T) — V(S,T) = c(S,T) = A(S). 


We observe that V is indeed a “down-and-out” barrier option. In other words, 1(down-and-in) plus 
1(down-and-out) equal to 1 vanilla call. This is because one and only one of the two barrier options 
can be active at expiry and whichever it is, its value is the value of a vanilla call. 


7.5 Asian options and lookback options 


In Asian options and lookback options, their payoff functions depend on the history of the underly- 
ing asset. For example, 


1. a European-type average strike option has the following payoff function 
1 rt 
max{Sp — = if S(r)dr, 0}. 
L Jo 
2. an American-type average strike option, 
1 t 
A(S,¢) = max{S — al S(r)dr, O}. 
0 
3. geometric mean, 
aDe 
A(S,T) = max{S — efo 8 5(7)4 gy, 
4. Lookback call, 
A(S,T) = max{S — J, 0}, LS Ae S(r) 


In general, the payoff depends on J, which is defined by 


r= f H8(0),7)¢r, 


where f is a smooth function. The payoff function is A(S, J). It is important to notice that I(t) is 
independent of S(t). The value of an asian option should depend on S, ¢ s well as I. Indeed, we 
shall see in the next chapter that df = fdt. The only randomness is through S, therefore V can be 
valued through a delta hedge. 

For the lookback option, it will be treated as a limiting case of an asian option. We shall discuss 
this in the next chapter. 
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Chapter 8 


Path-Dependent Options 


8.1 Introduction 


If the payoff depends on the history of the underlying asset price, such an option is called a path- 
dependent option. The Asian options and the Russian options (Lookback options) are the typical 
examples. The payoff functions for these options are, for example, 


1. average strike call option: A = max{S — 4 - S(r)dr, 0}, 
2. average rate call option: A = max{ + ie S(r)dr — E,0} 
3. geometric mean: the arithmetic mean z fi S(r)dr above is replaced by elo lox(S(r))ar, 


4. lookback strike put: A = max{maxo<;<7 S(T) — S, O}. 


5. lookback rate put: A = max{E — maxo<;<7 S(T), O}. 


8.2 General Method 


Let f be a smooth function, define 


me) = | PUSS Cr) ar 


In previous examples, f(.S(7),7) = S(r) for arithmetic mean and f(S(7),7) = log S(r) for 
geometric mean. 

Notice that J(t) is a random variable and is independent of S(t). (This is because J(t) is the sum 
of increment of functions of S' before time t, and each increment of S(T), 7 < t, is independent 
of S(t).) Therefore, we should introduce another independent variable I besides S to value the 
derivative V(S, J, t). 
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The stochastic differential equations governed by S' and J are 


° pdt + odz, 


dI(t) = f(S(t),t)dt. 


Notice that there is no noize term in dJ. The only randomness is from dS. Therefore, we can use 
delta hedge to eliminate this randomness. Namely, we consider the portfolio 


T=VAAS, 


as before. We have 


1 2 
dl = dV — AdS = (V, + 50°S? oat | VidI + VgdS — AdS. 


We choose A = a to eliminate the randomness in dII. From the arbitrage assumption, we arrive 


dene . GOV? _ GOV 
57 3 752 t fae r(V —-S—). 


i DS 


8.3 Average strike options 


8.3.1 European calls 


Let us consider an average strike call option with European exercise feature. Its payoff function is 
defined by 


1 i 
max{S — 7 i S(r)dr, 0). 


Or, in terms of I, A(S,I,T) = max{S — 4, O}. 
We notice that the modified Black-Scholes equation 


OV oo? 49 0°V OV 
Y%t+S bes t rg V =0, 
sae) ee hae ae 
and the initial data for the average strike options are invariant under the transformation: (S,I) > 
A(S, 1). Therefore, we expect that its solution is a function of the scale-invariant variable R = I/S. 


Notice that this is also reflected in that 
dR = (14+ (o? — p)R)dt — o Raz, 


depends on R only. Since the initial data can be expressed as A = S max{1 — R/T,0}, we may 
also expect that V = SH(R,t). This reduces one independent variable. 
Plug V = SH into the above modified Black-Scholes equation: 


OH OPH OH a) 


1 2 2 | | vd 
SH + 50°S' (255 + Say) + 5: SS +rS55(SH) — (SH) = 0. 
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From 
oO — OR dO _ RA 
OS ~=—s AS’ €@R-~~—SC SOR 
Oo? O10. . of? 0" 1 3) 5 oe 
aS ~ Stan’ stam ~ g2?Pap + ® op): 
we obtain 
1S 5 ROH. ¢ 1 OH | 90°H 
Sit + 50 S26 SOR? LS gzthap | DR2) 


+8? “Hp +rSH + rS7( ain rSH) =0 


Finally, we arrive 


Hy, + ii 1G +(1 
PS EO” OR? 
The payoff function 


A(R, T) = max{l1 - . 0} = H(R,T), (final condition) . 


we should require the boundary conditions. 


e H(oo,t) =0. Since as R > co implies S — 0, then V > 0, then H — 0. 


111 


e H(0,t) is finite. When R = 0, we have S(7) = 0, for all 7 with probability 1. That implies 


that A = 0, and consequently, H is finte at (0, t). 


Next, we expect that the solution is smooth up to R =. This implies that Hp, Hprp are finite at 
(0, ¢). We have the following two cases: (i) If R? aes O(1) 4 0as R — 0 then H = O(log R) 


ORz 


for R — 0. Or (ii) if Re = O(1) 4 0as R > O, then H = O(log R). Both cases contradict 
to H(0,t) being finite. Hence, we have RH p(0,t), R?HRR(0,t) are zeros as R — 0. Hence the 


boundary condition H (0, t) is finite is equivalent to H;(0,t) + Hp(0,t) = 0. 


This equation with boundary condition can be solved by using the hypergeometric functions. 


However, in practice, we solve it by numerical method. 


8.3.2 American call options 


We consider the average strike call option with American exercise feature. In this case, 


= 
5 DR? +(1-—rR) aR = 0 
H-A > 0 

1 9 520°H OH 2 
(A; 4 5° R DR? | ( rR) 5p )(H A) 0 Q, 


where A(R, t) = max{1— £,0}, R(t) = I(t)/S(t), I(t) = fo S(r)dr. 
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8.3.3. Put-call parity for average strike option 


We study the put-call parity for average strike options with European exercise feature. Consider a 
portfolio is C' — P. The corresponding payoff function is 


R R R 


Since the Black-Scholes equation in linear in R, we only need to solve the equation with final 
condition (i) H(R,T) = 1, (ii) H(R,T) = —2z. For (i), H(R,t) = 1. For (ii), since the final 
condition and the P.D.E. is linear in R, we expect that the solution is also linear in R. Thus, we 
consider #7 is of the following form 


A(t, R) = a(t) + d(t)R 


Plug this into the equation, we obtain 


and q q 
Het e=9 Gee 


with a(T’) = 0, b(T’) = —#. This differential equation can be solved easily: 
Le pied) I —r(T-t) 
b(t) = — Fe , a(t) =——(l-e HE 


Consequently, we obtain the put-call parity: 


1 1 Loaf 
= = = Gs g te 2p oS 
C—P S(1 apt e re s |, S(r)dr) 
t 
Ga ag ET ae ns i S(r)dr 
T Jo 


8.4 Lookback Option 


A lookback option is a derivate product whose payoff depends on the maximum or minimum of 
its underlying asset price. For instance, the payoff function for a lookback option with European 


exercise feature is 
VS — S(T). 
ae S(r) — S(T) 
Such an option is relatively expansive because it gives the holder an extremely advantageous payoff. 
As before, let us introduce J(t) = maxo<;<; S(T). Since S(7),7 < t are independent of S(t), 
we see that J(t) is independent of S(t). This suggest that we should introduce another independent 
variable J to value the lookback option in addition to S and t. We can derive a stochastic differential 
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equation for J as before. Indeed, it is dJ = 0. However, we shall give a more careful approach. We 
shall use the fact that for a continuous function S(-), 


‘lr 
t n 
guax |S) = lim. (f istr)rar)”. 


We leave its proof as an exercise. 
Remark. For the minimum, we have 


Let us introduce 


The s.d.e. for Jn, 


t+dt 4 t i 
a / (S(r))"dr)* — ( / (5(r))"dr)3 
= ee 


n 


Now, as before we consider the delta hedge: 
Il = P—AS. 
From the arbitrage assumption, we can derive the equation for P(S, Jn, t): 


1 S" OP 1 0?P 


Tl = Prdt4 dt + —0?S? 
d dt + aor yt + 50S Gantt 
OP 
= r(P——S)dt 
r( ng oa 


Taking n — oo, using the facts that J, + J and 5 < 1, we arrive 


Lg 0k OP - 
57 9 aap trae rP=0. 


Pot 


This is the usual Black-Scholes equation. The role of J here is only a parameter. This is consistent 
to the fact that 


dJ'= 0; 
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8.4.1. A lookback put with European exercise feature 
The range for Sis0 < S < J. This is because S < J, forO <¢ < T. We claim that 
P(0, J,t) = Je), 


Firstly, we have that A(0, J,T) = max{J — $,0} = J. Secondly, if S(t) = 0, then S(7r) = 
fort < + < T. The asset price process becomes deterministic. Therefore, the value of P is the 
discounted payoff: P(0, J,t) = Je". 


Next, we claim that 
OP 


OJ 
From js > 0, the current maximum cannot be the final maximum with probability 1. The value of 
P must be insensitive to a small change of J. 
We can use method of image to solve this problem. Its solution is given by 


S 


(J, J,t) = 


P= S(-1+N(dr)(1+8"1)) + eT ON (ds) — k-(5)*N (ds), 
where 
A o? 
ds = In(Q)-(r- SYP -dfovT = 
o2 
ds = [n(5)~(r- L(V 1)\/ovT 4 
o2 
dy = (nS) - (r+ SP -O|/ovT 
ia Bar 


8.4.2 Lookback put option with American exercise feature 
We have the following linear complementary equation, 
LpsP <0, P-A>0, (LesP)(P—- A) =0, 
where 
O° oS? 0" 6) 
= ters 
2 Os? Os 


Ts 


The final condition 
PUS J TV = NST) SS SS, 
The boundary condition 
OP 
Be By 
ap bat) = 
We require P, 2 oe ge 7 are continuous. 
For lookback call option, we simply replace maxo<;<;.5(7) by ming<;<;.5(7). For instance, 
its payoff is A = S(T) — mino<;<r S(T). 


Chapter 9 


Bonds and Interest Rate Derivatives 


9.1 Bond Models 


A bond is a long-term contract under which the issuer promises to pay the bondholder coupon 
payment (usually periodically) and principal (at the maturity dates). If there is no coupon payment, 
the bond is called a zero-coupon bond. The principal of a bond is called its face value. 


9.1.1 Deterministic bond model 


The value of a bond certainly depends on the interest rate. Let us first assume that the interest rate is 
deterministic temporarily, say r(T), t < 7 < T, is known. Let B(t, 7’) be the bond value at ¢ with 
maturity date T, k(t) be its coupon rate. This means that in a small dt, the holder receives coupon 
payment k(t) dt. From the no-arbitrage argument, 


dB + k(t) dt = r(t)Badt. 
together with the final condition: 

B(T,T) = F(face value), 
the bond value can be solved and has the following expression: 

T d T ae d 
B(t,T) =e fh "a LF + / k(r)els 74% ar] 
t 

Thus, the bond value is the sum of the present face value and the coupon stream. 


However, the life span of a bond is long (usually 10 years or longer), it is unrealistic to assume 
that the interest rate is deterministic. In the next subsection, we shall provide a stochastic model. 


9.1.2 Stochastic bond model 
Let us assume that the interest rate satisfies the following s.d.e.: 


dr = u(r, t)dt + w(r, t)dz, 
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where dz is the standard Wiener process. The drift wu and the variance w? are proposed by many 


researchers. We shall discuss these issues later. To find the equation for B with stochastic property 
of r, we consider a portfolio containing bonds with different maturity dates: 


Il = Vit T;) = AVE Yr, T2) = Vi cme AV). 
The change dII in a small time step dt is 


dil = dV, — AdVs, 


where 
o = pdt + o;dz, 
wi = 7 (Vie + wir + FwVior) 
= a wVie 


We we choose A = V,,,/V2,,, then the random term is canceled in dII. From the no-arbitrage 
argument, dII = rIdt. We obtain 


bLaiVi dt — ApeVedt = r(Vy — Av?) dt. 


This yields 
(Hi — 7)Vi/Vir = (M2 — 7) V2/Va,r, 


or equivalently 


PLE =, PO 
Ol 02 ; 
Since the left-hand side is a function of 7, while the right-hand side is a function of T2. Therefore, 
it is independent of T’. Let us express it as a known function X(r, t): 


uU-Tr 
= t). 
— = Xr.) 


Plug jz; and o; back to this equation, and drop the index 2, we obtain 
1 » 
Yt sw Ver + (u—Aw)V, —rV = 0. 


The function \(r,t) = 4 — is called the market price, since it gives the extra increase in expected 
instantaneous rate of return on a bond per an additional unit of risk. 

This stochastic bond model depends on three parameter functions u(r, t), w(r,t) and A(r,t). In 
the next section, we shall provide some model to determine them. 
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9.2 Interest models 


There are many interest rate models. We list some of them below. 

e Merton (1973): dr = adt + odz. 

e Vasicek (1977): dr = B(a — r)dt + odz. 

e Dothan (1978): dr = ordz. 

e Marsh-Rosenfeld (1983): dr = (ar?! + Br)dt + or?/*dz. 

e Cox-Ingersoll-Ross (1985) dr = B(a —r)dt + orfdz. 

e Ho-Lee (1986): dr = a(t)dt + odz. 

e Black-Karasinski (1991): dlnr = (a(t) — b(t) nr) dt + odz. 
The main requirements for an interest rate model are 

© positivity: r(t) > 0 almost surely, 

e mean reversion: 7 should tends to increase (or to decrease) and toward a mean. 
The C-I-R and B-K models have these properties. 

Below, we shall illustrate a unified approach proposed by Luo, Yen and Zhang. 
9.2.1 A functional approach for interest rate model 


The idea is to design r to be a function of x(t) and t, (ie. r = r(ax(t),t)) with x(t) governed by 
a simple stochastic process. We notice that the Ornstein-Uhlenbeck process dx = —nadt + adz 
has the property to tend to its mean (which is 0) time asymptotically. While the Bessel process 
dx = ¢«/xdt + odz has positive property. We then design the underlying basic process is the sum of 
these two processes: 


dx = (=n + =) dt + odz. 


In general, we allow 7, €, and o are given functions of t. With this simple process, we can choose 
r =r(a(t),t). Then all interest models mentioned above correspond to different choices of r(x, t), 
n, € and €. 


e Merton’s model: we choose € = 7 = 0,r= a+ at. 
e Dothan’s model: € = 7 =0,r = eto [2b 


e Ho-Lee: « =7=0,r=a2+ f5 a(s) ds. 


e Vasicek:€ =0,n=68,r=ax+a. 


C-LR: 8 = 2n, a = (0? + 2e)/(8n), r = ro 
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e Black-Karasinski: € = 0, r = exp(g(t, z)), where g = x + fo a(s) ds, 7 = U(t). 
With the interest rate model, the zero-coupon bond price V is given by 
VG@tyeoeE (ee penis) ae ae w) t<T. 


From the Feymann-Kac formula, V satisfies 


6) Lo €\ O 
ata : ( ua a Ox |v =o. 


This model depends three parameter functions e(t),7(t), a(t), and r(z,t). There is no unified 
theory available yet with this approach and the approach of the previous subsection. 


9.3 Convertible Bonds 


A convertible bond is a bond plus a call option under which the bond holder has the right to convert 
the bond into a common shares. Thus, it is a function of r, S,¢ and T’. Let the stochastic processes 
governed by S' and r are 


< = wpudt+oadzgs, 
dr = udt+wdz,. 


Suppose the correlation between dzg and dz p is 
dzg dz, = p(S,r, t)dt. 


The final value of the convertible bond V(r, S,T) = F’, the face value of the bond. Suppose the 
bond can be converted to nS at any time priori to 7. Then we have 


Vie Sat). > 8. 
We also have the boundary conditions: 


lin VS, by) Ss nS 


Soo 


lim: Vir, S,0)) =" 0; 


TCO 


At S = 0 orr = 0, we should require V(r, 0,t) or V(0, 5, t) to be finite. 
The Black-Scholes analysis for a convertible bond is similar to the analysis for a bond. Let 
V; = V(r, S,t, T;), i = 1, 2. Consider a portfolio 


T= AV, + AgV2 + AgS. 
In a small time step dt, the change of dII is 


dil = A;dV, + AgdV2 + AgsdS, 
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where ai 
— = pdt + op dz, + og dz, 
V; 
1 cme w 
He = (Viet 2S Vass + pSwVise + 2 Vier + poVis + Uae |, 
7 
1 
osi = FooMis 
7 
1 
Ori = 77 wir. 
7 
This implies 


dil = (AipaiVy ae Aope2V2 + Aus) dt 
+ (Ajog iV + Agogs2V2 + AoS) dzg 
+ (AioriVi + A2oy,2V2) dzp. 


We choose A,, Ay and As to cancel the randomness terms dz, and dzg. This means that 


(AjogiMi + Aoos2V2 + AoS) dzg = 0 
(Ajo;1V1 + Ao07,2V2) dz, = 0. 
And it yields 
dil = (AipiVy + Aopi2V2 + Aus) dt 
— r(AiVy + AoVe + AS) dt. 
Or equivalently, 


Ai(u1 — 7)Vy + Ao(pe2 — r)V2 + A(u— 7) S = 0. 


This equality together with the previous two give that there exist A, and Ag such that 


(1 = r) = Asosi le Ar Or] 
(2 = r) = Asos2 Tr Ar Or,2 
(u—Tr) = Aso 


The functions \, and Ag are called the market prices of risk with respect to r and S, respectively. 
Plug the formulae for ju;, 7; and os ;, we obtain the Black-Scholes equation for a convertible bond: 


2 2 2 2 92 
Oo A 0 sec Le eo 


a 
a 2> ase +’ asa + 2 Or gaa ee ae 


Or 
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Appendix A 


Basic theory of stochastic calculus 


A.1 From random walk to Brownian motion 


This section gives a more rigorous approach to construct Brownian motion from random walk. 
Without loss of generality, we shall construct motion B(t) for t € [0,1]. Let us consider a hierar- 
chical grids: for > 1 integer, define 


age 
A hierarchical independent random variables Af are Li.d.: At ~A, 


a 1 prob. 1/2 
~~ | -1- prob. 1/2. 
This corresponds to a binomial trial with sample space 


O :={(al, vos Oe) [45 =lor -1} 


The o-algebra F“ is defined to be F* := 2°" There are 2(2) elements in 9’. The Probability P* at 
each sample in 1° is defined to be Q- 2 
We define the random variables 


i< 

,=—) At. 7=1..., 2° 

vj: iz > J Per ess 
os 


and define a path x which is piecewise linear and its values at grid points are defined to be 


a! (t*) = a 
Each sample w € 0° corresponds to a unique piecewise linear path. We may also imagine the 
sample space 1 is the collection of such paths. 
The grids G’ := {t' a the path x’, the sample space 9 and the o-algebra F“ have the nested 
structure: 
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dG oe wiht =e. 


(2) Of! c Of with O54 = a; = a 
OVE er 
(4) The probability P’ on F*, as restricted to F‘~! is identical to P! 


With this nested structure, we can define the Brownian motion to be 
VSO. Fal 
é é 
with probability P := P* on F*. 
The Brownian motion B(t) has the following properties: 


1. B(O) = 0 and B(t) is continuous. Indeed, we should expect that B(-) is Hélder continuous 
continuous of order 1/2; 


2. If ty < te < ts < tg, then B(t,) — B(tz) and B(t2) — B(t1) are independent; 


3. B(s +t) — B(s) ~ N(0,t). This is mainly due to the central limit theorem. 


A.2. Brownian motion 

Let (Q, F, P) be a probability space. A process is a function X : [0,00) x (Q, F) + (£, 8B), such 
that for each t > 0, X(t) is a random variable. Here, E = R‘, B is the set of all Borel sets in E. 
Let 


Pe =O Missa) 
FF’ => of{X(s\33 >t} 


A process is called Markov if 
P(A| F,) = P(A| X(t), VA € F*. 
This is equivalent to 
P{X(r) € BJF} = P{X(r) € BIX(t)} Vr >t, 
A Markov process is characterized by its transition probability: 
P(t,z,s, B) := P{X(s)€ B| X(t) =z} 


with initial distribution 
P{X(0) € B} = v(B). 
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Theorem 1.10. Jf X is a Markov process, then the corresponding transition probability P satisfies 
Chapman-Kolmogorov equation: 


| Plto.to.t1, der) P(t, 21,t2,B) dx, = P(to, Xo, tz, B). 


Conversely, if P is a function satisfies Chapman-Kolmogorov equation, then there is a Markov 
process whose transition probability is P. 


Two standard Markov processes are the Wiener process and the Poisson process. The Wiener 
process has the transition probability density function 


lc —y|? 


= — 1)” 


exp( 


p(t, £, 8,y) = : 
(2r(t — s))4/2 
Such a distribution is called a normal distribution with mean x and variance s — t. 
Definition 1.8. Brownian motion: A process is called a Brownian motion(or a Wiener process) if 
1. B,— B, has normal distribution with mean 0 and variance t — s, 
2. it has independent increments: that is B, — Bs is independent of B, for allu < s < t, 
3. B, is continuous in t. 


It is easy to see that B, is markovian and its transition probability is 


1 exp( |x = yl? 
2n(t — s))4/2 s—t 


, 


p(t, £, 8, y) a ( 


Definition 1.9. A process { X;, | t > 0} is called martingale if 
(i) EX, < o, 
(ii) E(Xy|X.,0<s<th=X. 


This means that if we know the value of the process up to time t and X; = z, then the future 
expectation of X, is x. 


Theorem 1.11. I. Be is a martingale. 
2. B? —t is amartingale. 


Proof. 1. E(B,) = 0. From the fact that B;;, — B;, is independent of B;, for all s > 0, we 
obtain F'(B,4; — By | By) = 0, for all s > 0. Hence, 


E(Birs| Bi) = E((Birs — Be) + By | Bi) 
E((Bri+s — Br) | Be) + E(B: | Bz) 
Bs 
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2B (B) Src, 
3. Use 


Bi. = ((Bets — Bi) + Br)? 
= (Bis — Br)? +2Bi(Birs — Br) + B? 


and the fact that B,;, — B, is independent of B;, we obtain 
B(B},, | Br) = 8 + BP. 


Hence, 
E(B?,, — (t+) | By) = B?-t. 


We can show that the Brownian motion has infinite total variation in any interval. This means 
that 


lim)” |B(ti) — B(ti-1)| = 00. 
However, its quadratic variation, defined by 


[B, B](a, b) = lim)” |B(ti) — B(ti-a))? 


is finite. Here, a = to < ti <--: < tn = bisa partition of (a,b), and the limit is taken to be 
max(t; = 3) —> 0. 


Theorem 1.12. We have |B, B|(0,t) = t almost surely. 


Proof. Let us partition (0,t) evenly into 2” subintervals. Let T,, = an |B(t;) — B(ti-1)|?. We 
see that 


ET, = )~ E(\B(ts) — Bi-1))) = Soltis -— tial st. 


Var(In) = Var(S— |B(ti) — B(ti-1)|?) 
>= Var((B(ti) — B(i-1))”) 
2S 0 (ts — tin)? 


Dem 


Hence, \°°°_, Var(TI;,) < co. From Fubini theorem, 


E (Se, = en) < 00. 
n=1 


This implies E((T,, — ET;,)*) — 0 and hence, T,, + ET;, almost surely. 
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A.3 Stochastic integral 


We shall define the integral 
t 
[ tan, 
0 


The method can be applied with B replaced by a martingale, or a martingale plus a function with 
finite total variation. We shall require that f € H?, which means: 


(i) f(t) depends only on the history F; of B,, for s < t, 
(ii) [) El f/? < 00. 
For this kind of functions, we can define its It6’s integral as the follows. 


1. We define It6’s integral for f € H? and being step functions. Its It6’s integral is define by 
t n 
[ 1) 4B) =D MeayB) - BH). 
0 i=0 


2. We use above step functions f,, to approximate general function f € H?. Using the fact that, 
for step functions g € H?, 


E (( il @ ays) | = / “Bla(s)[2ds, 


one can show that 


almost surely. 


An typical example is 
: 1 t 
| OU C te Cue 
0 2 2 


From the definition, the integral can be approximated by 


In = > B(ti-1)(B(ti) — B(ti-a). 
i=l 


We have 
In = = [B? (ti) = B?(t;_1) — (B(ti) = B(ti-1))?] 
i=1 
= a) - hw) — Be? 
i=l 


We have seen that the second on the right-hand side tends to st almost surely. 
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A.4_ Stochastic differential equation 
A stochastic differential equation has the form 


A Markov process X is said to be a strong solution of this s.d.e. if it satisfies 
t t 
X,— Xp = | a(X5, 5) ds +f b(X 5, s)dB(s). 
0 0 
Theorem 1.13 (It6’s formula). Jf X satisfies the s.d.e. dX = adt + bdB, then 
1 
df (X(t) = (af'(X(t) + 56° f"(X(O))dt + F(X H)baB(t), 


We shall give a brief idea of the proof. In a small time step (t;_1, t;), let AB = B(t;)— B(ti_1), 
At = t; — t;_1. We have 


F(X (tia) + AX) — F(X (ta) = FAX + 5 f"(AX)? + o((AX)). 
We notice that 


(AX)? (aAt + BAB)? 


a” (At)? + 2abAt + b?(AB)? 
b?At + o(At). 


2 


Plug AX and (AX)? into the Taylor expansion formula for f. This yields the It6’s formula. 


A.5 Diffusion process 
For a s.d.e.(I.1p, we define the operator 7; on functions f € L?(IR) by 
Tif = Exz(f(X(s)),s > t. 
The operator 7; maps L?(R) to itself. It is a semigroup, that is 
De Le Ty its 


This can be proven from the Markovian property of X(-). Indeed, in terms of the transition proba- 
bility density function p(t, x, s, y), 


Tif = / amar) ay. 
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For a semigroup T;, we define its generator as 


Lf:= igs 
to0 t 


From It6’s formula, 
al 
Lf :=af' + 50 fF". 


This follows from It6’s formula and E(f, g(s) dB(s) = 0). 
With a fixed s > t and f, we define 


u(x,t) = (Tif) (a,t) = Exa(f(X(s)). 


Theorem 1.14. /f X satisfies the s.d.e. (7.1), then the associate u(x,t) := Ez2(f(X(s)), s >t 
satisfies the backward diffusion equation: 


uz + Lu = 0, 
and has final condition u(s,x) = f(z). 
Proof. First, we notice that 
u(x,t—h) = Eye n(f(X(s)) 


Ex t-nEx a e(f (X(s)) 
Ex t—n(u(X (), t)) 


This shifts final time from s to t. Now, we consider 


u(x,t — A) —u(a,t) nth (u( X(t), t) — u(X(t — h), t)) 
= aL, Lu(X(s),t) ds 
— Lu(x,t) 


as h — 0+. Next, u(x,s —h) = E,5-nf(X(s)), we have X(s) + x as h + 0+ almost surely. 
Hence u(z,s — h) > f(x) ash > 0+. 


Since u can be represent as 


u(t) = i p(a,t,,y) f(y) dy, 


we obtain that p satisfies 
pit+ Lop =0, 
and 
p(s, a, 8,y) = O(a — y). 
For diffusion equation with source term, its solution can be represented by the following Feymann- 
Kac formula. 
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Theorem 1.15 (Feymann-Kac). If X satisfies the s.d.e. (1.1), then the associate 


(et) i= Fas [x(s)) ex | g(X(r),T) ar| {sot (1.2) 
t 
solves the backward diffusion equation: 
ut + Lu+ gu= 0, 
with final condition u(s,x) = f(x). 


Proof. As before, we shift final time from s to f: 


s 


u(x,t — h) 


Exen|FX())exp [ g(X(0),7) dr 
= Bxp-nExee [£X()) exp fo X(7),7)ar| 


oe | Pye, | f(X(s)) exp i: Pecan ar| ee lexp / g(X(r),7) ar|] 


—h 
t 


+ Bas [uxt Nee | 


Pec as) ar| ; 


Here, we have used independence of X in the regions (t — h, t) and (t, s). Now, from It6’s formula, 
we have 


t 


u(x,t — h) — u(a, t) High ¥, luxe, t) exp | g(X(r),7) dr — u(X(t— h), 4] 


—h 


| 
S 
ie 
a 
—— 
Q 
—/~ 
= 
>< 
& 
ot 
oO 
ia 
Ke) 
—_ 
cS 
p= 
bas 
— 
= 
Q 
= 
nN—_” 


From this, it is easy to see that 


u(x, t- h) = u(a, t) = Lu(z, t) + g(x, Lule, ae 


